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Abstract — We  address  the  percolation-based  connectivity  of 
large-scale  ad  hoc  heterogeneous  wireless  networks,  where  sec¬ 
ondary  users  exploit  channels  temporarily  unused  by  primary 
users  and  the  existence  of  a  communication  link  between  two  sec¬ 
ondary  users  depends  on  not  only  the  distance  between  them  but 
also  the  transmitting  and  receiving  activities  of  nearby  primary 
users.  We  introduce  the  concept  of  connectivity  region  defined  as 
the  set  of  density  pairs — the  density  of  secondary  users  and  the 
density  of  primary  transmitters  —  under  which  the  secondary 
network  is  connected.  Using  theories  and  techniques  from  con¬ 
tinuum  percolation,  we  analytically  characterize  the  connectivity 
region  of  the  secondary  network  and  reveal  the  tradeoff  be¬ 
tween  proximity  (the  number  of  neighbors)  and  the  occurrence 
of  spectrum  opportunities.  Specifically,  we  establish  three  basic 
properties  of  the  connectivity  region — contiguity,  monotonicity  of 
the  boundary  and  uniqueness  of  the  infinite  connected  component, 
where  the  uniqueness  implies  the  occurrence  of  a  phase  transition 
phenomenon  in  terms  of  the  almost  sure  existence  of  either  zero 
or  one  infinite  connected  component;  we  identify  and  analyze 
two  critical  densities  which  jointly  specify  the  profile  as  well  as 
an  outer  bound  on  the  connectivity  region;  we  study  the  im¬ 
pacts  of  secondary  users’  transmission  power  on  the  connectivity 
region  and  the  conditional  average  degree  of  a  secondary  user 
and  demonstrate  that  matching  the  interference  ranges  of  the 
primary  and  the  secondary  networks  maximizes  the  tolerance  of 
the  secondary  network  to  the  primary  traffic  load.  Furthermore, 
we  establish  a  necessary  condition  and  a  sufficient  condition  for 
connectivity,  which  lead  to  an  outer  bound  and  an  inner  bound  on 
the  connectivity  region. 

Index  Terms — Cognitive  radio,  connectivity  region,  continuum 
percolation,  critical  densities,  heterogeneous  wireless  network, 
phase  transition. 


I.  Introduction 

THE  communication  infrastructure  is  becoming  increas¬ 
ingly  heterogeneous,  with  a  dynamic  composition  of  in¬ 
terdependent,  interactive  and  hierarchical  network  components 
with  different  priorities  and  service  requirements.  One  example 
is  the  cognitive  radio  technology  [1]  for  opportunistic  spectrum 
access  which  adopts  a  hierarchical  structure  for  resource  sharing 
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[2] .  Specifically,  a  secondary  network  is  overlaid  with  a  pri¬ 
mary  network,  where  secondary  users  identify  and  exploit  tem¬ 
porarily  and  locally  unused  channels  without  causing  unaccept¬ 
able  interference  to  primary  users  [2] . 

A.  Connectivity  and  Connectivity  Region 

While  the  connectivity  of  homogeneous  ad  hoc  networks  con¬ 
sisting  of  peer  users  has  been  well  studied  (see,  for  example, 

[3] — [10]),  little  is  known  about  the  connectivity  of  heteroge¬ 
neous  networks.  The  problem  is  fundamentally  different  from 
its  counterpart  in  homogeneous  networks.  In  particular,  the  con¬ 
nectivity  of  the  low-priority  network  component  depends  on  the 
characteristics  (traffic  pattern/load,  topology,  interference  toler¬ 
ance,  etc.)  of  the  high-priority  component,  thus  creating  a  much 
more  diverse  and  complex  design  space. 

Due  to  the  hierarchical  structure  of  spectrum  sharing,  a  com¬ 
munication  link  exists  between  two  secondary  users  if  the  fol¬ 
lowing  two  conditions  hold:  (Cl)  they  are  within  each  other’s 
transmission  range;  (C2)  they  see  a  spectrum  opportunity  de¬ 
termined  by  the  transmitting  and  receiving  activities  of  nearby 
primary  users  (see  Section  II-B-I).  We  say  that  a  topological  link 
exists  between  two  secondary  users  if  only  condition  (Cl)  holds. 

Using  theories  and  techniques  from  continuum  percolation, 
we  analytically  characterize  the  percolation-based  connectivity 
of  the  secondary  network  in  a  large-scale  ad  hoc  heterogeneous 
network.  Specifically,  we  consider  a  Poisson  distributed  sec¬ 
ondary  network  overlaid  with  a  Poisson  distributed  primary 
network  in  an  infinite  two-dimensional  Euclidean  space1.  We 
define  network  connectivity  as  the  existence  of  an  infinite 
connected  component  formed  by  communication  links  almost 
surely  (a.s.),  i.e.,  the  occurrence  of  percolation.  This  connec¬ 
tivity  characterizes  the  accessibility  of  two  secondary  users  at 
a  given  instant  and  is  thus  also  called  instantaneous  mutual 
connectivity.  Given  the  transmission  power  and  the  interfer¬ 
ence  tolerance  of  both  the  primary  and  the  secondary  users,  the 
connectivity  of  the  secondary  network  depends  on  the  density 
of  secondary  users  [due  to  (Cl)]  and  the  traffic  load  of  primary 
users  [due  to  (C2)]. 

We  thus  introduce  the  concept  of  connectivity  region  C,  de¬ 
fined  as  the  set  of  density  pairs  (As,  Xpr)  under  which  the  sec¬ 
ondary  network  is  connected,  where  As  denotes  the  density  of 
the  secondary  users  and  A pt  the  density  of  primary  transmitters 
(representing  the  traffic  load  of  the  primary  users).  As  illustrated 
in  Fig.  1,  a  secondary  network  with  a  density  pair  (As,  A  pt)  in¬ 
side  the  shaded  region  is  connected:  the  secondary  network  has 

'This  infinite  network  model  is  equivalent  in  distribution  to  the  limit  of  a 
sequence  of  finite  networks  with  a  fixed  density  as  the  area  of  the  network  in¬ 
creases  to  infinity,  i.e.,  the  so-called  extended  network  [11].  It  follows  from  the 
arguments  similar  to  the  ones  used  in  [12,  Chapter  3]  for  homogeneous  ad  hoc 
networks  that  this  infinite  ad  hoc  heterogeneous  network  model  represents  the 
limiting  behavior  of  large-scale  networks. 
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Fig.  1.  Connectivity  region  C  (the  upper  boundary  ApT(As)  is  defined  as  the 
supremum  density  of  the  primary  transmitters  that  ensures  connectivity  with  a 
fixed  density  of  the  secondary  users;  the  critical  density  \*s  of  the  secondary 
users  is  defined  as  the  infimum  density  of  the  secondary  users  that  ensures  con¬ 
nectivity  under  a  positive  density  of  the  primary  transmitters;  the  critical  density 
ApT  of  the  primary  transmitters  is  the  supremum  density  of  the  primary  trans¬ 
mitters  that  ensures  connectivity  with  a  finite  density  of  the  secondary  users). 

a  giant  connected  component  which  includes  infinite  number  of 
secondary  users.  The  existence  of  the  giant  connected  compo¬ 
nent  enables  bidirectional  communications  between  distant  sec¬ 
ondary  users  via  multihop  relaying.  On  the  other  hand,  a  sec¬ 
ondary  network  with  a  density  pair  (As,  Xpt)  outside  this  re¬ 
gion  is  not  connected:  the  network  is  separated  into  an  infinite 
number  of  finite  connected  components.  Consequently,  any  sec¬ 
ondary  user  can  only  communicate  with  users  within  a  limited 
range. 

It  is  shown  in  [13]  and  [14]  that  when  the  primary  network  is 
static,  this  instantaneous  connectivity  is  equivalent  to  the  con¬ 
nectivity  defined  in  terms  of  the  finiteness  of  the  minimum  mul¬ 
tihop  delay  (referred  to  as  fd-connectivity);  when  the  primary 
traffic  has  temporal  dynamics  (no  matter  how  small  the  range 
of  the  dynamics  is),  this  instantaneous  connectivity  is  stronger 
than  the  fd-connectivity  and  the  secondary  network  can  be  in¬ 
termittently  connected  as  long  as  it  is  topologically  connected 
[13],  [14]. 

The  objective  of  this  paper  is  to  establish  analytical  charac¬ 
terizations  of  the  connectivity  region  and  to  study  the  impact  of 
system  design  parameters  (in  particular,  the  transmission  power 
of  the  secondary  users)  on  network  connectivity.  The  main  re¬ 
sults  are  summarized  in  the  subsequent  two  subsections. 

B.  Analytical  Characterizations  of  the  Connectivity  Region 

We  first  establish  three  basic  properties  of  the  connectivity  re¬ 
gion:  contiguity,  monotonicity  of  the  boundary  and  uniqueness 
of  the  infinite  connected  component.  Specifically,  based  on  a 
coupling  argument,  we  show  that  the  connectivity  region  is  a 
contiguous  area  bounded  below  by  the  As -axis  and  bounded 
above  by  a  monotonically  increasing  function  ApT(As)  (see 
Fig.  1),  where  the  upper  boundary  ApT(As)  is  defined  as 

ApT(As)  =  sup{ApT  :  G(Xs,Xpp)  is  connected}  (1) 

with  G(Xs,Xpp)  denoting  the  secondary  network  of  density  As 
overlaid  with  a  primary  network  specified  by  the  density  Xpp 
of  the  primary  transmitters.  The  uniqueness  of  the  infinite  con¬ 
nected  component  is  established  based  on  ergodic  theory  and 
certain  combinatorial  results.  It  shows  that  once  the  secondary 
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network  is  connected,  it  contains  a  unique  infinite  connected 
component  a.s. 

Second,  we  identify  and  analyze  two  critical  parameters  of 
the  connectivity  region:  AJ  and  A *PT.  They  jointly  specify  the 
profile  as  well  as  an  outer  bound  on  the  connectivity  region.  Re¬ 
ferred  to  as  the  critical  density  of  the  secondary  users,  A  £  is  the 
infimum  density  of  the  secondary  users  that  ensures  connectivity 
under  a  positive  density  of  the  primary  transmitters 

X*s  =  inf  {As  :  3Xpp  >  0  s.t.  Q(X  s,  Xpp)  is  connected}. 

(2) 

We  show  that  A  J  equals  the  critical  density  Ac  of  a  homogeneous 
ad  hoc  network  (i.e.,  in  the  absence  of  primary  users),  which 
has  been  well  studied  [15].  This  result  shows  that  the  “takeoff’ 
point  in  the  connectivity  region  is  completely  determined  by  the 
effect  of  proximity — the  number  of  neighbors  (nodes  within  the 
transmission  range  of  a  secondary  user). 

Referred  to  as  the  critical  density  of  the  primary  transmitters, 
ApT  is  the  supremum  density  of  the  primary  transmitters  that 
ensures  the  connectivity  of  the  secondary  network  with  a  finite 
density  of  the  secondary  users: 

ApT  =  supjApp  :  3Xs  <  oc  s.t.  G(Xs,  Xpp)  is  connected}. 

_  (3) 

We  obtain  an  upper  bound  on  ApT,  which  is  shown  to  be  achiev¬ 
able  in  simulations.  More  importantly,  this  result  shows  that 
when  the  density  of  the  primary  transmitters  is  higher  than  the 
(finite)  value  given  by  this  upper  bound,  the  secondary  network 
cannot  be  connected  no  matter  how  dense  it  is.  This  parameter 
App  thus  characterizes  the  impact  of  opportunity  occurrence 
on  the  connectivity  of  the  secondary  network:  when  the  den¬ 
sity  of  the  primary  transmitters  is  beyond  a  certain  level,  there 
are  simply  not  enough  spectrum  opportunities  for  any  secondary 
network  to  be  connected. 

Since  a  precise  characterization  of  the  upper  boundary 
App  (As)  of  the  connectivity  region  is  intractable,  we  establish 
a  necessary  and  a  sufficient  condition  for  connectivity  to  pro¬ 
vide  an  outer  and  an  inner  bound  on  the  connectivity  region. 
The  necessary  condition  is  expressed  in  the  form  of  the  condi¬ 
tional  average  degree  of  a  secondary  user  and  is  derived  by  the 
construction  of  a  branching  process.  The  sufficient  condition 
is  obtained  by  the  discretization  of  the  continuum  percolation 
model  into  a  dependent  site  percolation  model. 

C.  Impact  of  Transmission  Power  on  Connectivity:  Proximity 
Versus  Opportunity 

The  study  of  the  impact  of  the  secondary  users’  transmis¬ 
sion  power  ptx  on  network  connectivity  reveals  an  interesting 
tradeoff  between  proximity  and  opportunity  in  the  design  of 
heterogeneous  networks.  As  illustrated  in  Fig.  2,  we  show 
that  increasing  ptx  enlarges  the  connectivity  region  C  along 
the  As -axis  (i.e.,  better  proximity  leads  to  a  smaller  “takeoff’ 
point),  but  at  the  price  of  reducing  C  along  the  Apt -axis. 
Specifically,  with  a  large  ptx,  few  secondary  users  experience 
spectrum  opportunities  due  to  their  large  interference  footprint 
with  respect  to  the  primary  users.  This  leads  to  a  poor  tolerance 
to  the  primary  traffic  load  parameterized  by  Xpp. 

The  transmission  power  ptx  of  the  secondary  network  should 
thus  be  chosen  according  to  the  operating  point  of  the  hetero- 


REN  et  al.:  CONNECTIVITY  OF  HETEROGENEOUS  WIRELESS  NETWORKS 


4317 


Fig.  2.  Simulated  connectivity  regions  for  two  different  transmission  powers 
(ptx  denotes  the  transmission  power  of  the  secondary  users  and  the  large  ptx  is 
3"  times  the  small  ptx ,  where  a  is  the  path-loss  exponent). 

geneous  network  given  by  the  density  of  the  secondary  users 
and  the  traffic  load  of  the  co-existing  primary  users.  Using  the 
tolerance  to  the  primary  traffic  load  as  the  performance  mea¬ 
sure,  we  show  that  the  interference  range  rj  of  the  secondary 
users  should  be  equal  to  the  interference  range  Rp  of  the  pri¬ 
mary  users  in  order  to  maximize  the  upper  bound  on  the  critical 
density  A pT  of  the  primary  transmitters.  Given  the  interference 
tolerance  of  the  primary  and  secondary  users,  we  can  then  de¬ 
sign  the  optimal  transmission  power  ptx  of  the  secondary  users 
based  on  that  of  the  primary  users. 

D.  Related  Work 

To  our  best  knowledge,  the  connectivity  of  large-scale  ad  hoc 
heterogeneous  networks  has  not  been  characterized  analytically 
or  experimentally  in  the  literature.  There  are  a  number  of  classic 
results  on  the  connectivity  of  homogeneous  ad  hoc  networks. 
For  example,  it  has  been  shown  that  to  ensure  either  1 -connec¬ 
tivity  (there  exists  a  path  between  any  pair  of  nodes)  [3],  [4] 
or  ^-connectivity  (there  exist  at  least  k  node-disjoint  paths  be¬ 
tween  any  pair  of  nodes)  [5],  the  average  number  of  neighbors 
of  each  node  must  increase  with  the  network  size.  On  the  other 
hand,  to  maintain  a  weaker  connectivity — p-connectivity  (i.e., 
the  probability  that  any  pair  of  nodes  is  connected  is  at  least  p), 
the  average  number  of  neighbors  is  only  required  to  be  above  a 
certain  “magic  number’  which  does  not  depend  on  the  network 
size  [6]. 

The  theory  of  continuum  percolation  has  been  used  by 
Dousse  et  al.  in  analyzing  the  connectivity  of  a  homogeneous 
ad  hoc  network  under  the  worst-case  mutual  interference  [7], 
[8].  In  [9],  [10],  the  connectivity  and  the  transmission  delay  in 
a  homogeneous  ad  hoc  network  with  static  or  dynamic  on-off 
links  are  investigated  from  a  percolation-based  perspective. 

In  [13]  and  [14],  the  impact  of  temporal  dynamics  of  primary 
traffic  on  the  connectivity  and  multihop  delay  of  the  secondary 
network  is  studied.  In  particular,  it  is  shown  that  in  the  presence 
of  primary  traffic  dynamics,  the  secondary  network  is  connected 


(either  instantaneously  or  intermittently)  as  long  as  it  is  topolog¬ 
ically  connected.  Optimal  power  control  in  heterogeneous  net¬ 
works  has  been  studied  in  [16],  which  focuses  on  a  single  pair  of 
secondary  users  in  a  Poisson  network  of  primary  users.  The  im¬ 
pacts  of  secondary  users’  transmission  power  on  the  occurrence 
of  spectrum  opportunities  and  the  reliability  of  opportunity  de¬ 
tection  are  analytically  characterized. 

E.  Organization  and  Notations 

The  rest  of  this  paper  is  organized  as  follows.  Section  II 
presents  the  Poisson  model  of  the  heterogeneous  network.  In 
particular,  the  conditions  for  the  existence  of  a  communication 
link  in  the  secondary  network  is  specified  based  on  a  rigorous 
definition  of  spectrum  opportunity.  In  Section  III,  we  introduce 
the  concept  of  connectivity  region  and  establish  its  three  basic 
properties.  The  two  critical  densities  are  analyzed,  followed 
by  a  necessary  and  a  sufficient  condition  for  connectivity.  In 
Section  IV,  we  demonstrate  the  tradeoff  between  proximity 
and  opportunity  by  studying  the  impact  of  the  secondary  users’ 
transmission  power  on  the  connectivity  region  and  on  the 
conditional  degree  of  a  secondary  user.  Further,  we  obtain  the 
optimal  transmission  power  of  the  secondary  users,  adopting 
the  secondary  network’s  tolerance  of  the  primary  traffic  load  as 
the  performance  measure.  Section  V  contains  detailed  proofs 
of  the  main  results  and  Section  VI  concludes  the  paper. 

Throughout  the  paper,  we  use  capital  letters  for  parameters 
of  the  primary  users  and  lowercase  letters  for  those  of  the  sec¬ 
ondary  users. 

II.  Network  Model 

We  consider  a  Poisson  distributed  secondary  network  over¬ 
laid  with  a  Poisson  distributed  primary  network  in  an  infinite 
two-dimensional  Euclidean  space.  The  models  of  the  primary 
and  secondary  networks  are  specified  in  the  following  two  sub¬ 
sections. 

A.  The  Primary  Network 

The  primary  transmitters  are  distributed  according  to  a  two- 
dimensional  Poisson  point  process  with  density  A pt-  Each  pri¬ 
mary  transmitter’s  receiver  is  uniformly  distributed  within  its 
transmission  range  Rp.  Here  we  have  assumed  that  all  primary 
transmitters  use  the  same  transmission  power  and  the  trans¬ 
mitted  signals  undergo  isotropic  path  loss.  Based  on  the  dis¬ 
placement  theorem  [17,  Chapter  5],  it  is  easy  to  see  that  the  pri¬ 
mary  receivers  form  a  two-dimensional  Poisson  point  process 
with  density  Apt-  Note  that  the  two  Poisson  processes  formed 
by  the  primary  transmitters  and  receivers  are  correlated. 

B.  The  Secondary  Network 

The  secondary  users  are  distributed  according  to  a  two-di¬ 
mensional  Poisson  point  process  with  density  A 5,  independent 
of  the  Poisson  processes  of  the  primary  transmitters  and  re¬ 
ceivers.  The  transmission  range  of  the  secondary  users  is  de¬ 
noted  by  rp. 

Communication  Links:  In  contrast  to  the  case  in  a  homoge¬ 
neous  network,  the  existence  of  a  communication  link  between 
two  secondary  users  depends  on  not  only  the  distance  between 
them  but  also  the  availability  of  the  communication  channel 
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Fig.  3.  Definition  of  spectrum  opportunity. 

(i.e.,  the  presence  of  a  spectrum  opportunity).  The  latter  is  deter¬ 
mined  by  the  transmitting  and  receiving  activities  in  the  primary 
network  as  described  below. 

As  illustrated  in  Fig.  3,  there  exists  an  opportunity  from  A, 
the  secondary  transmitter,  to  B ,  the  secondary  receiver,  if  the 
transmission  from  A  does  not  interfere  with  nearby  primary  re¬ 
ceivers  in  the  solid  circle  and  the  reception  at  B  is  not  affected 
by  nearby  primary  transmitters  in  the  dashed  circle  [18].  Re¬ 
ferred  to  as  the  interference  range  of  the  secondary  users,  the 
radius  rj  of  the  solid  circle  at  A  depends  on  the  transmission 
power  of  A  and  the  interference  tolerance  of  the  primary  re¬ 
ceivers,  whereas  the  radius  Rj  of  the  dashed  circle  (the  interfer¬ 
ence  range  of  the  primary  users2)  depends  on  the  transmission 
power  of  the  primary  users  and  the  interference  tolerance  of  B. 

It  is  clear  from  the  above  discussion  that  spectrum  opportuni¬ 
ties  depend  on  both  transmitting  and  receiving  activities  of  the 
primary  users.  Furthermore,  spectrum  opportunities  are  a, sym¬ 
metric.  Specifically,  a  channel  that  is  an  opportunity  when  A 
is  the  transmitter  and  B  the  receiver  may  not  be  an  opportu¬ 
nity  when  B  is  the  transmitter  and  A  the  receiver.  In  other 
words,  there  exist  unidirectional  communication  links  in  the 
secondary  network.  Since  unidirectional  links  are  difficult  to  uti¬ 
lize  in  wireless  networks  [19],  we  only  consider  bidirectional 
links  in  the  secondary  network  when  we  define  connectivity.  As 
a  consequence,  when  we  determine  whether  there  exists  a  com¬ 
munication  link  between  two  secondary  users,  we  need  to  check 
the  existence  of  spectrum  opportunities  in  both  directions. 

To  summarize,  under  the  disk  signal  propagation  and  inter¬ 
ference  model,  there  is  a  (bidirectional)  link  between  A  and  B 
if  and  only  if  (Cl)  the  distance  between  A  and  B  is  at  most  rp\ 
(C2)  there  exists  a  bidirectional  spectrum  opportunity  between 
A  and  B,  i.e.,  there  are  no  primary  transmitters  within  distance 
Rj  of  either  A  or  B  and  no  primary  receivers  within  distance  rj 
of  either  A  or  B. 

Connectivity:  We  interpret  the  connectivity  of  the  secondary 
network  in  the  percolation  sense:  the  secondary  network  is  con¬ 
nected  if  there  exists  an  infinite  connected  component  a.s. 

Based  on  the  above  conditions  (Cl,  C2)  for  the  existence  of  a 
communication  link,  we  can  obtain  an  undirected  random  graph 
G(^s^pt)  corresponding  to  the  secondary  network,  which 
is  determined  by  three  Poisson  point  processes:  the  secondary 

2Since  the  minimum  transmission  power  for  successful  reception  is,  in  gen¬ 
eral,  higher  than  the  maximum  allowable  interference  power,  it  follows  that  the 
transmission  range  Rp  of  primary  users  is  smaller  than  R  j .  Furthermore,  under 
the  disk  signal  propagation  and  interference  model,  we  have  Rp  —  3RT  (0  < 
(3  <  1).  A  similar  relationship  holds  for  rp  and  rT. 
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Fig.  4.  Realization  of  the  heterogeneous  network.  The  random  graph 
Q(XsApt)  consists  of  all  the  secondary  nodes  and  all  the  bidirectional  links 
denoted  by  solid  lines.  The  solid  circles  with  radii  Ri  denote  the  interference 
regions  of  the  primary  transmitters  within  which  secondary  users  cannot 
successfully  receive  and  the  dashed  circles  with  radii  r7  denote  the  required 
protection  regions  for  the  primary  receivers  within  which  the  secondary  users 
should  refrain  from  transmitting. 

users  with  density  A 5,  the  primary  transmitters  with  density 
A pt  and  the  primary  receivers  with  density  A pt  (correlated 
with  the  process  of  the  primary  transmitters)3.  See  Fig.  4  for  an 
illustration  of  G(\s,  A  pt). 

The  question  we  aim  to  answer  in  this  paper  is:  when  is  the 
secondary  network  connected,  i.e.,  when  does  percolation  occur 
in  f/(As,  Apt)? 

III.  Analytical  Characterizations  of  the 
Connectivity  Region 

Given  the  transmission  power  and  the  interference  tolerance 
of  both  the  primary  and  the  secondary  users  (i.e.,  Rp ,  Rj ,  rp 
and  rj  are  fixed),  the  connectivity  of  the  secondary  network  is 
determined  by  the  density  As  of  the  secondary  users  and  the 
density  Apt  of  the  primary  transmitters.  We  thus  introduce  the 
concept  of  connectivity  region  C  of  a  secondary  network,  which 
is  defined  as  the  set  of  density  pairs  (As,  A  pt)  under  which  the 
secondary  network  <?(As,  Apt)  is  connected  (see  Fig.  1) 

C  =  {(As,  Apt)  •  f/(As,  Apt)  is  connected.}. 

A.  Basic  Properties  of  the  Connectivity  Region 

We  next  establish  three  basic  properties  of  the  connectivity 
region. 

Theorem  1:  Basic  Properties  of  the  Connectivity  Region. 
Tl.l  The  connectivity  region  C  is  contiguous,  that  is,  for 
any  two  points  (Asi,  Apti),  (As2?  A pt2 )  £  C,  there  exists 
a  continuous  path  in  C  connecting  the  two  points. 

T1.2  The  lower  boundary  of  the  connectivity  region  C  is 
the  As-axis.  Let  ApT(As)  denote  the  upper  boundary  of 
the  connectivity  region  C  defined  by  (1),  then  we  have  that 
A pT(Xs)  is  monotonically  increasing  with  As- 

3The  two  Poisson  point  processes  of  the  primary  transmitters  and  receivers 
are  essentially  a  snap  shot  of  the  realizations  of  the  primary  transmitters  and 
receivers.  In  different  time  slots,  different  sets  of  primary  users  become  active 
transmitters/receivers.  Thus,  even  if  a  secondary  user  is  isolated  at  one  time  due 
to  the  absence  of  spectrum  opportunities,  it  may  experience  an  opportunity  at  a 
different  time  and  be  connected  to  other  secondary  users. 
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(a)  (b) 


Fig.  5.  Two  realizations  of  the  Poisson  heterogeneous  network  (black  stars  denote  primary  transmitters,  green  plus  signs  denote  primary  receivers,  red  dots  denote 
secondary  users  and  blue  segments  denote  the  bidirectional  links  between  secondary  users).  We  have  removed  secondary  users  who  do  not  see  opportunities  for 
clarity.  The  simulation  parameters  are  given  by  Rp  —  50  m,  =  80  m,  Xs  —  650  km-2,  rp  =  50  m,  rT  —  80  m  and  the  critical  density  in  this  case  is 

Ac(50)  ft;  576  km-2  (a)  percolation  occurs  (APT  =  10  km-2)  (b)  percolation  does  not  occur  (APT  =  20  km-2). 


T1 .3  There  exists  either  zero  or  one  infinite  connected  com¬ 
ponent  in  <5(As,  Apt)  a.s. 

Proof:  The  proofs  of  Tl.l  and  T1.2  are  based  on  the 
coupling  argument,  a  technique  frequently  used  in  continuum 
percolation  [15,  Section  2.2].  The  proof  of  T  1.3  is  based  on 
the  ergodicity  of  the  random  model  driven  by  the  three  Poisson 
point  processes  of  the  primary  transmitters,  the  primary  re¬ 
ceivers  and  the  secondary  users.  The  details  of  the  proofs  are 
given  in  Section  V-B.  ■ 

Tl.l  and  T1.2  specify  the  basic  structure  of  the  connectivity 
region,  as  illustrated  in  Fig.  1.  T1.3  implies  the  occurrence  of  a 
phase  transition  phenomenon,  that  is,  there  exists  either  a  unique 
infinite  connected  component  a.s.  or  no  infinite  connected  com¬ 
ponent  a.s.  This  uniqueness  of  the  infinite  connected  component 
excludes  the  undesirable  possibility  of  having  more  than  one 
(maybe  infinite)  infinite  connected  component  in  the  secondary 
network.  We  point  out  that  such  a  property  is  not  always  present 
in  wireless  networks.  Two  examples  where  more  than  one  infi¬ 
nite  connected  component  exists  in  a  homogeneous  ad  hoc  net¬ 
work  can  be  found  in  [20] . 


B.  Critical  Densities 

In  this  subsection,  we  study  the  critical  density  A  J  of  the  sec¬ 
ondary  users  and  the  critical  density  XPT  of  the  primary  trans¬ 
mitters  defined  by  (2)  and  (3),  respectively.  We  have  the  fol¬ 
lowing  theorem. 


Theorem  2:  Critical  Densities. 

Given  Rp ,  Rj ,  rp  and  77,  we  have: 

T2.1  AJ  =  Xc(rp),  where  Xc(rp)  is  the  critical  density  for 
a  homogeneous  ad  hoc  network  with  transmission  range  rp 
(i.e.,  in  the  absence  of  the  primary  network). 

T2.2 


Ac(l) 


^ PT  4ma x{Rprj} 


(4) 


where  the  constant  Ac(l)  is  the  critical  density  for  a  homo¬ 
geneous  ad  hoc  network  with  a  unit  transmission  range. 


Proof:  The  basic  idea  of  the  proof  of  T2. 1  is  to  approximate 
the  secondary  network  Q(Xs ,  A pt)  by  a  discrete  edge-percola¬ 
tion  model  on  the  grid.  This  discretization  technique  is  often 
used  to  convert  a  continuum  percolation  model  to  a  discrete 
site/edge  percolation  model  (see,  for  example,  [15,  Chapter  3], 
[8]).  The  details  of  the  proof  are  given  in  Section  V-C-I. 

The  proof  of  T2.2  is  based  on  the  argument  that  if  there  is  an 
infinite  connected  component  in  the  secondary  network,  then  an 
infinite  vacant  component  must  exist  in  the  two  Poisson  Boolean 
models  driven  by  the  primary  transmitters  and  the  primary  re¬ 
ceivers,  respectively.  The  key  point  is  to  carefully  choose  the 
radii  of  the  two  Poisson  Boolean  models  in  order  to  obtain  a 
valid  upper  bound  on  ApT.  The  details  of  the  proof  can  be  found 
in  Section  V-C-II.  ■ 

Fig.  5(a)  shows  one  realization  of  the  heterogeneous  network 
when  Xs  is  slightly  larger  than  Xc(rp )  and  A  pt  is  small.  At  least 
one  left-to-right  (L-R)  crossing  and  at  least  one  top-to-bottom 
(T-B)  crossing  can  be  found  in  the  square  network.  It  is  thus 
expected  that  these  L-R  and  T-B  crossings  in  finite  square  re¬ 
gions  can  form  an  infinite  connected  component  in  the  whole 
network  on  IR2.  If  we  slightly  increase  A  pt,  then  we  observe 
from  Fig.  5(b)  that  the  reduction  in  spectrum  opportunities  elim¬ 
inates  numerous  communication  links  in  the  secondary  network, 
creating  several  disjoint  small  components. 

Fig.  6  shows  a  simulation  example  of  the  connectivity  region, 
where  the  upper  bound  on  the  critical  density  X*PT  of  the  pri¬ 
mary  transmitters  given  in  T2.2  appears  to  be  achievable. 


C.  A  Necessary  Condition  for  Connectivity 

In  this  subsection,  we  establish  a  necessary  condition  for  con¬ 
nectivity  which  is  given  in  terms  of  the  average  conditional  de¬ 
gree  of  a  secondary  user.  This  condition  agrees  with  our  intu¬ 
ition:  the  secondary  network  cannot  be  connected  if  the  degree 
of  every  secondary  user  is  small. 

Let  I  (A,  d,  rx/tx)  denote  the  event  that  there  exist  primary  re¬ 
ceivers/transmitters  within  distance  d  of  a  secondary  user  A.  Let 
I(v4,  d ,  rx/tx)  denote  the  complement  of  I  (A,  d ,  rx/tx).  Since  a 
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Fig.  6.  Simulated  connectivity  regions  when  rp  —  150  m,  rj  =  240  m, 
Rp  —  1 00  m  and  RT  —  120  m.  The  blue  dashed  line  is  the  upper  bound  (4)  on 
the  critical  density  A *PT  of  primary  transmitters  given  in  T2.2.  The  area  of  the 
simulated  heterogeneous  network  is  2000  m  x  2000  m.  For  a  fixed  density  As 
of  the  secondary  users,  the  upper  boundary  AJ>,r(A<?)  is  equal  to  the  minimum 
density  of  the  primary  transmitters  such  that  over  all  the  1000  realizations,  the 
percentage  of  the  ones  in  which  there  exists  at  least  one  L-R  crossing  is  below 
50%.  The  intuitive  reason  for  choosing  the  existence  of  an  L-R  crossing  as  the 
criterion  for  connectivity  is  illustrated  in  Fig.  5  and  is  discussed  in  Section  V-C. 

secondary  user  is  isolated  if  it  does  not  see  a  spectrum  opportu¬ 
nity,  we  focus  on  secondary  users  who  experience  spectrum  op¬ 
portunities  and  define  the  conditional  average  degree  r  of  such 
a  secondary  user  A  as 

r  =  E[deg(A)\  l(A,r/,rx)  fl  I  (A, i2j,tx)]  (5) 

where  deg(A)  denotes  the  degree  of  A,  rj  the  interference  range 
of  the  secondary  users  and  Rj  the  interference  range  of  the  pri¬ 
mary  users.  Notice  that  the  degree  of  A  is  the  number  of  sec¬ 
ondary  users  within  the  transmission  range  of  A  and  experi¬ 
encing  bidirectional  opportunities  with  A.  We  arrive  at  the  fol¬ 
lowing  necessary  condition  for  connectivity. 

Theorem  3:  A  necessary  condition  for  the  connectivity  of 
G(Xs,  Apt)  is  r  >  1,  where  r  is  the  conditional  average  degree 
of  a  secondary  user  defined  in  (5). 

Proof:  The  basic  idea  is  to  construct  a  branching  process, 
where  the  conditional  average  degree  r  is  the  average  number  of 
offspring.  This  branching  process  provides  an  upper  bound  on 
the  number  of  secondary  users  in  a  connected  component.  If  r  < 
1,  then  the  branching  process  is  finite  a.s.  It  thus  follows  that 
there  is  no  infinite  connected  component  a.s.  in  <?(As,Apt)- 
Details  can  be  found  in  Section  V-D.  ■ 

To  apply  the  necessary  condition  given  in  Theorem  3,  the 
conditional  average  degree  r  of  a  secondary  user  A  needs  to 
be  evaluated  based  on  the  network  parameters.  Let  B  be  a  sec¬ 
ondary  user  randomly  and  uniformly  distributed  within  trans¬ 
mission  range  rp  of  A.  Let  g(\pT,rp,ri,  Rp,  Ri)  denote  the 
probability  of  a  bidirectional  opportunity  between  A  and  B  con¬ 
ditioned  on  the  event  that  A  sees  an  opportunity.  Based  on  the 
statistical  equivalence  and  independence  of  different  points  in 
a  Poisson  point  process,  the  conditional  average  degree  r  of  a 
secondary  user  A  is  given  by  this  conditional  probability  g(-)  of 
a  bidirectional  opportunity  between  A  and  a  randomly  chosen 


neighbor  multiplied  by  the  average  number  of  neighbors  of  A, 
i.e., 

H  =  (Xsirrl)  ■  g(XpT,  rp,  r/,  Rp ,  i?/).  (6) 

The  detailed  derivation  of  (6)  and  an  expression  for  gf)  are 
given  in  Appendix  A.  It  is  also  shown  in  Appendix  A  that  g(-)  is 
a  strictly  decreasing  function  of  Apt-  Thus  g_1(-),  the  inverse 
of  g(-)  with  respect  to  Apt,  is  well-defined. 

Combining  (6)  with  Theorem  3,  we  obtain  an  outer  bound  on 
the  connectivity  region.  Specifically,  let  r(Xs ,  Apt)  denote  the 
conditional  average  degree  of  a  secondary  user  in  G(\s,  Apt)- 
Then  those  density  pairs  (As,  Xpp)  satisfying  //(As,  A pT)  <  1 
are  outside  the  connectivity  region. 

Corollary  1:  Given  Rp,  Rj,  rp  and  rj9  an  outer  bound  on  the 
connectivity  region  C  is  given  by 

xpt = 9-1  ( A§) 

where  g_1(-)  is  the  inverse  of  the  conditional  probability  g(-) 
with  respect  to  Apt- 

D.  A  Sufficient  Condition  for  Connectivity 

In  this  subsection,  we  establish  a  sufficient  condition  for  con¬ 
nectivity,  which  provides  an  inner  bound  on  the  connectivity  re¬ 
gion  and  a  criterion  for  checking  whether  a  secondary  network 
is  connected. 

Our  sufficient  condition  for  connectivity  is  established  by 
using  the  discretization  technique.  The  continuum  percolation 
model  is  mapped  onto  a  dependent  site-percolation  model  C 
in  the  following  way.  As  illustrated  in  Fig.  7,  we  partition  IR2 
into  (dashed)  squares  with  side  length  ds  and  locate  a  site  at  the 
center  of  each  square.  Sites  whose  associated  dashed  squares 
share  at  least  one  common  point  are  considered  connected  (as 
illustrated  by  solid  lines  in  Fig.  7).  Thus,  each  site  is  connected 
to  eight  neighbors4  (see  the  eight  neighbors  0\ , . . . ,  Og  of  site 
O  in  Fig.  7).  Let  Bo  be  the  associated  dashed  square  of  O,  then 
O  is  occupied  if  there  exists  in  Bo  at  least  one  secondary  user 
who  sees  an  opportunity. 

Since  the  largest  distance  between  two  points  in  two  neigh¬ 
boring  dashed  squares  is  2xf2ds,  it  follows  that  if  we  set  ds  = 
then  every  pair  of  secondary  users  in  two  neighboring 
dashed  squares  are  within  transmission  range  rp  of  each  other. 
Based  on  the  definitions  of  occupied  site  in  C  and  communi¬ 
cation  link  in  the  secondary  network,  we  conclude  that  the  ex¬ 
istence  of  an  infinite  occupied  component  (a  connected  compo¬ 
nent  consisting  of  only  occupied  sites)  in  C  implies  the  existence 
of  an  infinite  connected  component  in  the  secondary  network. 

Due  to  the  fact  that  spectrum  opportunities  are  spatially  de¬ 
pendent,  the  state  of  one  site  is  correlated  with  the  states  of  its 
adjacent  sites.  Thus,  the  above  site-percolation  model  C  is  a  de¬ 
pendent  model.  Define  the  dependence  range  k  as  the  minimum 
distance  such  that  the  state  of  any  two  sites  at  distance  larger 
than  k  are  independent,  where  the  distance  between  two  sites 

4For  the  commonly  used  square  site-percolation  model,  each  site  has  four 
neighbors.  The  site-percolation  model  constructed  here  provides  a  better  inner 
bound. 
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Fig.  7.  Illustration  of  the  dependent  site-percolation  model  C  with  side  length 
cl s  (solid  dots  denote  sites,  solid  lines  denote  edges  connecting  every  two  sites 
and  dashed  lines  denote  the  squared  partition). 

is  the  minimum  number  of  neighboring  sites  that  must  be  tra¬ 
versed  from  one  site  to  the  other.  Then  the  dependence  range  of 
C  is  given  by 


Corollary  2:  A  sufficient  condition  for  the  connectivity  of 
<?(As,  A pT)  is 


A  pt  < 


q\(2Hl)2 

v37 


where  I(Ri ,  Rp,  rj)  is  defined  in  (8)  and  k  is  the  dependence 
range  of  the  site-percolation  model  defined  in  (7). 


IV.  Impact  of  Transmission  Power:  Proximity 
Versus  Opportunity 

In  this  section,  we  study  the  impact  of  the  secondary  users’ 
transmission  power  on  the  connectivity  and  the  conditional  av¬ 
erage  degree  of  the  secondary  network.  As  has  been  illustrated 
in  Fig.  2,  there  exists  a  tradeoff  between  proximity  and  oppor¬ 
tunity  in  designing  the  secondary  users’  transmission  power  for 
connectivity.  Specifically,  increasing  the  transmission  power  of 
the  secondary  users  leads  to  a  smaller  critical  density  AJ  of  the 
secondary  users,  but  at  the  same  time,  a  lower  tolerance  to  the 
primary  traffic  load  manifested  by  a  smaller  critical  density  \*PT 
of  the  primary  transmitters. 

A.  Impact  on  the  Conditional  Average  Degree 


k  =  8max|i?/  +  ^,77  +  /rp~|  -  1.  (7) 

Let  pc  denote  the  upper  critical  probability  of  C  which  is  de¬ 
fined  as  the  minimum  probability  of  occupancy  p*  such  that  if 
the  probability  of  occupancy  p  >  p* ,  an  infinite  occupied  com¬ 
ponent  containing  the  origin  exists  in  C  with  a  positive  prob¬ 
ability  (wpp.).  Since  the  dependence  range  k  of  C  is  finite,  it 
follows  from  [12,  Theorem  2.3.1]  that  pc  <  1.  Now  we  present 
a  sufficient  condition  for  connectivity  in  the  following  theorem. 

Theorem  4:  Let  pc  denote  the  upper  critical  probability  of  the 
dependent  site-percolation  model  C  specified  above.  Define 

I(r,  Rp,  rj)  =  2  jf  (8) 

where  Sj(t,  Rp,rj)  is  the  common  area  of  two  circles  with  radii 
Rp  and  rj  and  centers  t  apart.  Then  the  secondary  network  is 
connected  if  the  equation  shown  at  the  bottom  of  the  page  holds. 


As  discussed  in  Section  III-C,  the  expression  for  the  condi¬ 
tional  average  degree  p  can  be  decomposed  into  the  product 
of  two  terms:  As7rr^  and  g(\pT,  RP,  Ri)-  The  first 
term  is  the  average  number  of  neighbors  of  a  secondary 
user,  which  increases  with  the  transmission  power  ptx  of  the 
secondary  users  (i.e.,  enhanced  proximity).  The  other  term 
^(Apt,  rp,ri,Rp,  Rj)  is  the  conditional  probability  of  a  bidi¬ 
rectional  opportunity,  which  decreases  with  ptx  due  to  reduced 
spectrum  opportunities.  This  tension  between  proximity  and 
opportunity  is  illustrated  in  Fig.  8,  where  we  observe  that 
the  impact  of  ptx  on  proximity  dominates  when  ptx  is  small 
(p  increases  with  ptx)  while  its  impact  on  the  occurrence  of 
opportunities  dominates  when  ptx  is  large  ( p  decreases  with 
Ptx)' 

Corollary  3:  Let  ptx  be  the  transmission  power  of  secondary 
users  and  p  the  conditional  average  degree  defined  in  (5),  then 
under  the  disk  signal  propagation  and  interference  model  we 
have5 


Proof:  The  proof  is  based  on  the  ergodicity  of  the  hetero¬ 
geneous  network  model  and  its  relation  with  the  constructed 
dependent  site-percolation  model  C.  Details  can  be  found  in 
Section  V-E.  ■ 

By  applying  a  general  upper  bound  on  the  upper  critical  prob¬ 
ability  pc  for  a  site-percolation  model  with  finite  dependence 
range  [12,  Theorem  2.3.1],  we  arrive  at  the  following  corollary. 


M  =  O  (( ptx )  as  ptx  — >  oo 

where  a  is  the  path-loss  exponent. 

Proof:  We  show  this  corollary  by  deriving  an  upper  bound 
on  p.  Details  can  be  found  in  Appendix  B.  ■ 

5Here  we  use  the  Big  O  notation:  f(  x  )  =  0(g(x))  as  x  — ►  oo  if  and  only  if 
3  M  >  0,  ■V o  >  0  such  that  |/(.t)|  <  M\g(x)\  for  all  x  >  x0. 


1  —  exp 


8 


>  pc. 
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Fig.  8.  Conditional  average  degree  jl  of  secondary  users  vs  transmission  range 
rp  of  secondary  users.  rp  oc  (ptx)  “  »  where  ptx  is  the  transmission  power  of 
secondary  users  and  a  is  the  path-loss  exponent  and  simulation  parameters  are 
given  by  XPT  —  2.5  km-2,  Rp  —  200  m,  Rj  —  250  m,  A^  =  25  km-2, 
n  —  rp/ 0.8. 

For  a  homogeneous  network,  the  average  degree  of  a  user  is 
\7TVp,  which  increases  with  ptx  at  rate  (ptx)2^a-  In  sharp  con¬ 
trast,  this  corollary  tells  us  that  for  a  heterogeneous  network, 
when  ptx  is  large  enough,  the  conditional  average  degree  p  of 
a  secondary  user  actually  decreases  with  ptx  at  least  as  fast  as 

G Ptx)~2/a ■ 

B.  Impact  on  the  Connectivity  Region 

From  the  scaling  relation  of  the  critical  density  [15,  Propo¬ 
sition  2.11],  we  know  that  in  a  homogeneous  two-dimensional 
network 


—  2  —  — 

^■c(Xp)  —  Ac(l)  (r p)  oc  ( Ptx )  a 

where  the  constant  Ac(l)  is  the  critical  density  for  a  homoge¬ 
neous  ad  hoc  network  with  a  unit  transmission  range.  Thus, 
if  each  secondary  user  adopts  a  high  transmission  power,  then 
Xc(rp)  decreases.  It  follows  from  T2.1  that  the  critical  den¬ 
sity  A *s  of  secondary  users  required  to  achieve  connectivity  de¬ 
creases  due  to  the  enhanced  proximity. 

On  the  other  hand,  from  the  upper  bound  on  the  critical  den¬ 
sity  A pT  of  the  primary  transmitters  given  in  T2.2,  we  have  that 

X*PT  =  O  ((ptx)-2/“)  as  ptx  ->  oo 


Fig.  9.  Example  of  the  upper  bound  on  A *PT  as  a  function  of  rj  (parameters 
are  given  by  RT  —  120  m,  rp  —  0.625rj). 


leads  to  the  “largest”  connectivity  region  (largest  in  the  sense 
that  its  connectivity  region  contains  all  regions  achievable  with 
any  finite  transmission  power  ptx  of  the  secondary  users).  Thus, 
the  optimal  design  of  ptx  depends  on  the  operating  point  of  the 
heterogeneous  network.  For  instance,  when  a  sparse  secondary 
network  is  overlaid  with  a  primary  network  with  low  traffic  load, 
a  large  ptx  may  be  desirable  to  achieve  connectivity.  The  oppo¬ 
site  holds  when  a  dense  secondary  network  is  overlaid  with  a 
primary  network  with  high  traffic  load. 

Focusing  on  a  secondary  network  whose  density  exceeds  the 
critical  density  of  the  secondary  users,  we  address  the  design 
of  its  transmission  power  for  the  maximum  tolerance  to  the  pri¬ 
mary  traffic.  Due  to  its  tractability  and  achievability,  the  upper 
bound  on  the  critical  density  A PT  of  the  primary  transmitters 
given  in  T2.2  is  used  as  the  performance  measure  (see  Fig.  6  for 
a  simulation  example). 

Theorem  5:  Let  rj  and  Rj  denote  the  interference  range  of 
the  secondary  and  the  primary  users.  For  a  fixed  Rj ,  the  upper 
bound  on  XPT  given  in  T2.2  is  maximized  when  the  primary  and 
secondary  networks  have  matching  interference  ranges:  rj  = 
Rt- 

Proof:  Since  under  the  disk  signal  propagation  and  inter¬ 
ference  model,  rp  =  [3rj  for  some  [3  G  (0, 1),  the  upper  bound 
on  A pT  can  be  written  as 


A 


* 

PT 


< 


Xc  ( 1 ) 

4R]-{32rj 

(1) 

(4—/32)r2 


for  n  <  Ri 
for  rj  >  Rj. 


where  we  have  assumed  that  rp  =  f3rj  for  some  (3  G  (0, 1) 
under  the  disk  signal  propagation  and  interference  model.  Thus, 
when  the  transmission  power  ptx  of  the  secondary  network  is 
large  enough,  the  critical  density  XPT  of  the  primary  transmit¬ 
ters  decreases  with  ptx  at  least  as  fast  as  (ptx)~2^a  due  to  re¬ 
duced  spectrum  opportunities. 

C.  Optimal  Design  of  Transmission  Power 

Due  to  the  tension  between  proximity  and  opportunity,  there 
does  not  exist  a  transmission  power  of  the  secondary  users  that 


Then  the  above  theorem  can  be  readily  shown  by  finding  the 
maximal  point  for  the  two  cases:  rj  <  Rj  and  rj  >  Rj.  ■ 

An  example  of  the  upper  bound  on  XPT  is  plotted  as  a  func¬ 
tion  of  rj  in  Fig.  9.  Notice  that  there  is  a  distinct  difference 
in  the  slope  on  the  two  sides  of  the  optimal  point.  As  a  conse¬ 
quence,  the  operating  region  of  rj  <  Ri  is  preferred  over  that  of 
r j  >  Rj  when  the  optimal  point  rj  =  Rj  cannot  be  achieved. 
We  point  out  that  the  desired  operating  region  of  rj  <  Rj  is  the 
typical  case  of  a  secondary  network  coexisting  with  a  privileged 
primary  network. 


REN  et  al. :  CONNECTIVITY  OF  HETEROGENEOUS  WIRELESS  NETWORKS 

V.  Proofs 

In  this  section,  we  present  proofs  of  the  main  results  presented 
in  Sections  III-IV.  We  start  with  a  brief  overview  of  several 
basic  results  in  percolation  theory  that  will  be  used  in  the  proofs. 

A.  Percolation  Theory:  Some  Basic  Results 

Poisson  Boolean  Model:  Poisson  Boolean  model  [15],  often 
referred  to  as  B(X,  p,  A),  is  specified  by  two  elements:  a 
Poisson  point  process  X  on  Ud  with  density  A  and  a  radius  p, 
a  random  variable  with  a  given  distribution.  Under  this  model, 
each  point  in  X  is  the  center  of  a  ball  in  Rd  with  a  random 
radius  distributed  according  to  the  distribution  of  p.  Radii  asso¬ 
ciated  with  different  points  are  independent  and  they  are  also 
independent  of  points  in  X.  Under  a  Poisson  Boolean  model, 
the  whole  space  is  partitioned  into  two  regions:  the  occupied 
region,  which  is  the  region  covered  by  at  least  one  ball  and  the 
vacant  region,  which  is  the  complement  of  the  occupied  region. 
We  define  occupied  (vacant)  components  as  those  connected 
components  in  the  occupied  (vacant)  region. 

Assume  that  nodes  in  a  homogeneous  ad  hoc  network  form 
a  Poisson  point  process  with  density  A  and  their  transmission 
range  is  r.  It  is  easy  to  see  that  the  connectivity  of  this  network 
can  be  studied  through  examination  of  the  occupied  connected 
components  in  the  corresponding  Poisson  Boolean  model 
B(X,  r/2 ,  A). 

Sharp  Transition  in  Two  Dimensions:  For  the  Poisson 
Boolean  model  in  two  dimensions,  the  sharp  transition  phe¬ 
nomenon  is  remarkable  in  the  sense  that  the  critical  density 
for  the  a.s.  existence  of  infinite  occupied  components  is  equal 
to  that  for  the  a.s.  existence  of  infinite  vacant  components. 
Let  Ac(2 p)  denote  the  critical  density  for  the  Poisson  Boolean 
model  B(X,  p,  A),  then  we  have  that: 

•  when  A  <  Ac(2p),  there  is  no  infinite  occupied  component 
a.s.  and  there  is  a  unique  infinite  vacant  component  a.s.; 

•  when  A  >  Ac(2p),  there  is  a  unique  infinite  occupied  com¬ 
ponent  a.s.  and  there  is  no  infinite  vacant  component  a.s. 

The  exact  value  of  Ac  is  not  known.  For  a  deterministic  radius  p, 
simulation  results  [21]  indicate  that  Ac(2p)  ^  0.36p-2,  while 
rigorous  bounds  0.192p-2  <  Ac(2p)  <  0.843p-2  are  provided 
in  [15]  and  [22]. 

Crossing  Probabilities:  A  continuous  curve  in  the  occupied 
region  is  called  an  occupied  path.  An  occupied  path  7  is  an  oc¬ 
cupied  L-R  crossing  of  the  rectangle  {0  <  x  <  li}  x  {0  < 
V  <  h}  if  7  intersects  both  the  left  and  the  right  boundaries 
of  the  rectangle,  i.e.,  7  n  {{x  =  0}  x  {0  <  y  <  l2 })  ^  <t>, 
7  n  ({ x  =  h}  x  {0  <  y  <  I2 })  /  4>  and  the  segment  between 
the  two  intersecting  points  is  fully  contained  in  the  rectangle. 
Similarly,  we  define  an  occupied  T-B  crossing  by  requiring  that 
7  intersects  with  the  top  and  bottom  boundaries  of  the  rectangle. 
Let 

<r((h,h),  A,  L-R) 

=  Pr{3  an  occupied  L  —  R  crossing  of  [0,  l\ ]  x  [M2]} 
<r((h,h),  A,  T-B) 

=  Pr{3  an  occupied  T  —  B  crossing  of  [0,  li\  x  M} 
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>  ApT2  Ap 
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path  connecting  the  two  points 

(Api,  Appi) 

and  (A 52 7  ^pt2)  in 

the  connectivity 

region  C.  (a) 

Appi 

<  A  PT2  • 

(b)  A pti  >  A pt2- 

denote  the  two  crossing  probabilities  in  the  rectangle  [0,  Zi]  x 
[0,  Z2].  Then  for  a  Poisson  Boolean  model  B(X,  p,  A)  in  two 
dimensions  with  a.s.  bounded  p,  we  have  [15,  Corollary  4.1] 
that  for  any  k  >  1 

iil!mM„((U,B),A,L-R)  =  {j;  (9) 

Due  to  the  symmetry  of  the  Poisson  Boolean  model, 
similar  results  hold  for  the  T-B  crossing  probability 

cr((n,  kn ),  A,  T  —  B). 

B.  Proof  of  Theorem  1 

Proof  of  Tl.l:  To  prove  Tl.l,  it  suffices  to  show  that  for  any 
two  given  points  (Asi,  Xpti)  and  (As2,  App2)  i nC,  we  can  find 
a  path  in  C  that  connects  these  two  points.  In  particular,  the  path 
we  constructed  is  given  by  a  horizontal  segment  and  a  vertical 
segment  as  shown  in  Fig.  10,  where  we  assume,  without  loss  of 
generality,  that  Asi  <  As2- 

Consider  case  (a)  in  Fig.  10  where  Appi  <  App2 •  Case 
(b)  can  be  proven  similarly.  First,  we  show  every  point 
(As,  Xpti)  (Asi  <  As  <  As2)  on  the  horizontal  segment  be¬ 
longs  to  C.  Let  A'  =  As  —  Asi.  A  Poisson  point  process  X  with 
density  As  is  statistically  equivalent  to  the  superposition  of  a 
Poisson  point  process  X\  with  density  Asi  and  an  independent 
Poisson  point  process  X'  with  density  A'.  It  follows  that  any 
realization  of  the  heterogeneous  network  with  densities  As  and 
Xpti  can  be  generated  by  adding  more  secondary  nodes  to  a 
realization  of  the  heterogeneous  network  with  densities  Asi  and 
ApTi-  Thus,  the  existence  of  an  infinite  connected  component 
in  5(Asi,  Appi)  implies  the  existence  of  an  infinite  connected 
component  in  C/(As,  Appi).  We  thus  have  that  (As,  Appi)  £  C 
for  (Asi  <  As  <  As2)- 

Now  we  know  that  the  two  end  points  (As2?Appi)  and 
(As2?  A pp2)  of  the  vertical  segment  belong  to  C.  For  a  point 
(As2?  App)  (Appi  <  A pp  <  App2)  on  the  vertical  segment, 
let  A'  =  A pt2  ~  App,  then  any  realization  of  the  heteroge¬ 
neous  network  with  densities  As2  and  App  can  be  obtained 
by  independently  removing  each  primary  transmitter-receiver 
pair  with  probability  X' /Xpp2  from  a  realization  of  the  het¬ 
erogeneous  network  with  densities  As2  and  A pp2.  It  follows 
from  the  definition  of  communication  link  in  the  secondary 
network  (see  Section  II-B-I)  that  the  existence  of  an  infinite 
connected  component  in  <?(As25App2)  implies  the  existence 
of  an  infinite  connected  component  in  <?(As,App).  Thus,  we 
have  (As2?App)  £  C  (Appi  <  App  <  App2). 
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Proof  of  Theorem  1.2:  Suppose  that  (As,  App)  £  C  (Xpp  > 
0) ,  then  by  using  the  coupling  argument  for  showing  that  the  ver¬ 
tical  segment  belongs  to  C  in  the  above  proof  of  T  1.1,  we  con¬ 
clude  that  (As,  0)  E  C ,  i.e.,  the  As-axis  is  the  lower  boundary 
of  C. 

Suppose  that  As2  >  Xsi  >  0.  In  order  to  prove  the  mono¬ 
tonicity  of  ApT(As)  with  As  it  suffices  to  show  that  V  Xpp  >  0, 
if  (Asi,  Apt)  £=  C  then  (As2,  Xpp)  E  C.  This  is  a  direct  conse¬ 
quence  of  the  coupling  argument  for  showing  that  the  horizontal 
segment  belongs  to  C  in  the  above  proof  of  T  1.1. 

Proof  of  Theorem  1.3 :  We  first  establish  the  ergodicity6  of 
the  heterogeneous  network  model. 

Lemma  1:  The  heterogeneous  network  model  is  ergodic. 

Proof  of  Lemma  1:  The  proof  of  this  lemma  is  inspired 
by  the  proof  of  the  ergodicity  of  the  Poisson  Boolean  model 
[15,  Proposition  2.8].  The  difficulty  here  is  that  for  the  hetero¬ 
geneous  network  model,  we  have  two  correlated  Poisson  point 
processes:  the  primary  transmitters  and  the  primary  receivers. 
The  definition  of  the  shift  transformation  for  the  primary  net¬ 
work  model  is  thus  more  complicated  than  the  standard  Poisson 
Boolean  model  with  a  deterministic  radius  p.  To  prove  Lemma 
1,  we  first  show  the  ergodicity  of  the  primary  network  model 
and  then  we  show  the  mixing  property7  of  the  secondary  net¬ 
work  model.  Since  the  primary  network  model  is  independent 
of  the  secondary  network  model,  it  follows  from  [24,  Theorem 
2.6.1]  that  the  heterogeneous  network  model  is  ergodic. 

Let  X  denote  the  Borel  cr-algebra  in  Rd  and  TV  the  set  of  all 

d 

simple  counting  measures8  on  X  .  Construct  a  cr-algebra 
for  TV  generated  by  sets  of  the  form 

{n  E  TV  :  n(A )  =  k} 
d 

where  A  E  X  and  k  is  an  integer.  A  point  process  X  can  now 
be  defined  as  a  measurable  mapping  from  a  probability  space 
(12,  X  ,  P)  into  (TV,  pV)  [23,  Chapter  7].  The  measure  p  on 
pL  induced  by  X  is  defined  as  p(G)  =  P(X~1(G)),  for  all 
G  E  c/K. 

In  order  to  define  the  shift  transformation  on  12,  it  is  conve¬ 
nient  to  identify  (12,  X  )  with  (TV,  ,,  L).  Let  u(A)  denote  the 

d 

number  of  points  ini  E  J  ,  V  a;  E  O  and  Tx  be  the  shift 
according  to  a  vector  x  E  Ud.  Then  Tx  induces  a  shift  transfor¬ 
mation  Sx  :  O  — ^  12  through  the  following  equation  for  every 

Ae  md: 

(Sxlj)(A)  =  cop-1  A).  (10) 

6  A  measure  preserving  (m.p.)  dynamical  system  (fl,  X  ,  fi,  T)  consists  of 
a  set  Q,  a  o -algebra  A  of  measurable  subsets  of  Q,  a  nonnegative  measure 
H  on  (S2,  X  )  and  an  invertible  m.p.  transformation  T  :  O  — >  O  such  that 
ji(T~1F)  —  fi(F)  for  all  F  G  X  .  A  set  F  G  X  is  said  to  be  T-invariant  if 

T~XF  —  F.  The  m.p.  dynamical  system  (fl,  X  ,  fi,  T)  is  said  to  be  ergodic 
if  the  cr-algebra  of  T-invariant  sets  is  trivial,  i.e.,  for  any  invariant  set,  either  it 
has  measure  0  or  its  complement  has  measure  0. 

7An  m.p.  dynamical  system  (fl,  X  ,  T)  is  said  to  be  mixing  if  for  all 
E,F  e  X  ,  lim  n{TnE  n  F)  =  (i(E)fi(F). 

n  — >  oo 

d 

8A  simple  counting  measure  on  X  is  an  integer- valued  measure  for  which 
the  measures  of  bounded  Borel  sets  are  all  finite  and  the  measure  of  a  point  is 
at  most  1 . 


Let  (12 pt,  X  pp,  Ppt )  be  the  probability  space  of 
the  Poisson  point  process  Xpp  for  the  primary  transmit¬ 
ters  with  density  Xpp.  Let  12 pp  be  the  product  space 
rineN  rW  CRP  for  the  primary  receivers,  where 
Grp  =  {(x,y)  :  x2  +  y2  <  Rp}.  Then  we  equip  Opp  with 
the  usual  product  cr-algebra  and  product  measure  Ppr  with  all 
marginal  probability  measure  being  pu,  where  pu  is  a  uniform 
probability  measure  on  Crp  .  Finally,  we  set  12 p  =  12 pp  x  LLpr 
and  equip  12  p  with  the  product  measure  Pp  =  Ppp  x  Ppr  and 
the  usual  product  cr-algebra.  It  follows  that  the  primary  network 
model  is  a  measurable  mapping  from  12 p  into  Npp  x  LLpr 
defined  by  (ujpp,  upr)  — >(Xpp((jpp),  ljpr),  where  Npp  is 
specified  in  the  definition  of  the  point  process. 

The  positions  of  the  primary  transmitters  corresponding  to 
(wpt,  ojpr)  E  12 pp  x  LtpR  are  easily  known  from  ujpp. 
For  the  primary  receivers,  the  positions  are  obtained  as  follows. 
Consider  binary  cubes 
2 

K(n ,  z)  JJ(^2_n,  (zi  +  l)2_n]  for  all  n  E  N  and  z  E  Z2. 
i= 1 

For  each  primary  transmitter  xpp ,  there  exists  a  unique  smallest 
integer  no  =  no(xpp)  such  that  it  is  contained  in  a  binary  cube 
iT(no,  z (no,  xpp))  which  contains  no  other  primary  transmit¬ 
ters.  The  relative  position  of  xpp  ’  s  receiver  with  respect  to  xpp 
is  then  given  by  ujpr  (n0 , 2  (n0 ,  xPT ) ) . 

Let  ei,  e2  denote  the  unit  vectors  in  IR2,  then  the  translation 
Te.  :  [R2  — >  [R2  (z  =  1,  2)  defined  by  x  — >  x  +  ei  induces  a  shift 
transformation  Uei  on  12 pr  through  the  equation 

(■ UeiupR)(n ,  z)  =  ujpr(u,  z  -  2 71  ef),  for  i  =  1,2. 

Hence  Te.  also  induces  a  shift  transformation  Te.  on  Lip  = 
12 pp  x  LIpr  as  follows: 

Te>P)  =  ( Seiupp ,  UeiupR ),  fori  =  1,2 

where  Sej  is  defined  in  (10).  By  using  techniques  similar  to  the 
proof  of  Boolean  models  [15,  Proposition  2.8],  we  have  that  the 
m.p.  dynamical  system  (!2p,  X  p,  Pp,  Tei)  is  ergodic. 

Since  the  transmission  range  rp  of  secondary  users  is  fixed, 
the  probability  space  of  the  secondary  network  model  is  the 
probability  space  (12 5  ,&s,  Ps)  for  the  Poisson  point  process 
Xs  of  secondary  users  with  density  A5.  It  follows  from  the 
proof  of  Poisson  point  processes  [15,  Proposition  2.6]  that  the 
m.p.  dynamical  system  (!2s,  X  5,  Ps,  Sei)  is  mixing. 

Since  the  primary  network  model  is  independent  of  the  sec¬ 
ondary  network  model,  the  sample  space  of  the  heterogeneous 
network  model  12  can  be  written  as  the  product  of  12 p  and  Os, 
i.e.,  O  =  Op  x  Os-  We  equip  O  with  product  measure  P  = 
Pp  x  Ps  and  the  usual  product  cr-algebra.  Similarly,  the  trans¬ 
lation  Te.  (i  =  1,2)  induces  a  transformation  Te%  on  O  = 
Op  x  Os,  which  is  given  by 

Te.  (id)  =  (feiLUp,Se.LUs)- 

Then  it  follows  from  [24,  Theorem  2.6.1]  that  the  product  m.p. 
dynamical  system  (O,  X  ,  P,  Tei)  is  ergodic.  Since  the  cr-al¬ 
gebra  invariant  under  the  transformation  group  {Tz  :  z  E  I2} 
is  a  subset  of  the  cr- algebra  invariant  under  the  transformation 
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Tei,  we  conclude  that  {T2  :  z  G  Z2}  acts  ergodically,  i.e.,  the 
heterogeneous  network  model  is  ergodic.  ■ 

Based  on  Lemma  1,  we  have  the  following  lemma. 

Lemma  2:  The  number  of  infinite  connected  component  in 
Q(X 5,  A pt)  is  a  constant  a.s.  and  it  can  only  take  value  from 
(0,  1,  oo}. 

Proof  of  Lemma  2:  Let  K  denote  the  (random)  number  of 
infinite  connected  components  in  (/(As,  A  pt),  then  since  for  all 
k  >  0,  the  event  {K  =  k}  is  invariant  under  the  group  of  shift 
transformations,  it  follows  from  Lemma  1  and  the  definition  of 
ergodicity  that  Pr {K  =  k}  =  0  or  1.  Consequently,  we  have 
that  K  is  a.s.  constant.  Then  it  suffices  to  exclude  the  possibility 
of  K  >  2.  This  is  shown  by  contradiction,  that  is,  if  there  exist 
K  >  2  infinite  connected  components,  then  they  can  be  linked 
together  as  one  connected  component  wpp.  The  proof  is  inspired 
by  the  proof  of  Proposition  3.3  in  [15]  and  a  major  difference  is 
that  here  we  need  to  consider  the  impact  of  the  primary  network 
on  the  connectivity  of  the  secondary  network. 

Suppose  that  there  are  K  >  2  infinite  connected  components 
a.s.  If  we  remove  all  the  secondary  nodes  centered  inside  a  box 
B  =  [— n,n]2,  then  the  resulting  secondary  network  should 
contain  at  least  K  unbounded  components  a.s.  Let,  for  A  C  R2, 
G[A\  denote  the  graph  formed  by  secondary  nodes  in  A.  Given 
a  box  B  and  e  >  0,  consider  the  event  E(B,e)  shown  at  the 
bottom  of  the  page. 

Partition  the  box  B  into  squares  with  side  length  a  >  0 
and  let  <Sa  =  {Si, . . . ,  Sn}  denote  the  collection  of  all  the 
squares  which  are  adjacent  to  the  boundary  of  B.  Clearly,  for 
a  box  B  and  e  >  0,  we  can  find  a  =  a(B,e)  G  (0,rp/\/5) 
and  r]  =  77(a)  >  0  such  that  for  any  point  x  B  with 
d{x,B)  <  rp  —  e/2,  there  exists  a  square  S  =  S(x)  G  Sa 
for  which  we  have  supyGS  d(x,  y)  <  rp  —  rj.  This  means  that, 
if  we  center  in  each  square  of  Sa  a  secondary  node  and  there 
are  neither  primary  transmitters  nor  primary  receivers  within  a 
bigger  box  5  =  [— n  — max{r/,i?/},n  +  max{rj,i?/}]2,  then 
every  infinite  component  U  in  Q[BC ]  with  d(U,  B)  <  rp  —  e  is 
connected  to  some  secondary  node  in  Sa . 

Let  E(a,rj)  be  the  event  that  each  square  in  Sa  contains  at 
least  one  secondary  node  and  E(B)  the  event  that  there  are  nei¬ 
ther  primary  transmitters  nor  primary  receivers  within  B.  Since 
E{a,rj)  depends  on  the  configuration  of  secondary  nodes  inside 
the  box  B ,  E(B ,  e)  depends  on  the  configuration  of  secondary 
nodes  outside  B  and  the  configuration  of  primary  nodes,  based 
on  the  independence  of  the  primary  network  and  the  secondary 
network,  we  have 


Pr(E(B,e)  D  E(a,r])  C\  E(B)) 

=  Pr(E(B}  e))Pr (E(a,  t]))Py(E(B)\E(B,  e)). 
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If  E(B ,  e),  E(a,rj)  and  E(B)  all  occur  wpp.,  then  there  is 
only  one  infinite  connected  component9  wpp.  By  using  argu¬ 
ments  similar  to  the  proof  for  [15,  Proposition  3.3],  we  have 
that  there  exists  a  large  enough  box  B  and  e  >  0  such  that 
Pr (E(B,e))  >  0.  Obviously,  Pr(i?(a,  77))  >  0.  Moreover,  it 
is  easy  to  see  that  P{E{B)\E{B ,  e))  >  P{E{B))  >0.  ■ 

Now  we  have  that  the  number  K  of  infinite  connected  com¬ 
ponents  is  equal  to  zero,  one  or  infinity  a.s.  To  exclude  the  pos¬ 
sibility  of  K  =  00,  we  can  use  the  proof  of  [15,  Theorem  3.6] 
for  Poisson  Boolean  models,  which  is  based  on  several  combi¬ 
natorial  results.  The  details  are  omitted. 

C.  Proof  of  Theorem  2 

Proof  of  T2.1:  To  prove  T2.1,  it  suffices  to  show  that: 

(a)  for  any  A 5  <  A c(rp),  the  secondary  network  is  not  con¬ 
nected  for  any  A  pt  >  0; 

(b)  for  any  A 5  >  A c(rp),  there  exists  a  ApT(As)  >  0 
such  that  V  A  pt  <  ApT(As),  the  secondary  network  is 
connected. 

From  Section  IV-A,  we  know  that  for  a  Poisson  homogeneous 
ad  hoc  network  with  density  A  and  transmission  range  r ,  the  nec¬ 
essary  and  sufficient  condition  for  connectivity  is  A  >  Ac(r). 
Since  the  existence  of  an  infinite  connected  component  in  the 
secondary  network  implies  the  existence  of  an  infinite  connected 
component  in  the  homogeneous  ad  hoc  network  with  the  same 
density  and  the  same  transmission  range,  by  using  a  coupling 
argument,  we  conclude  that  when  A 5  <  A c(rp),  there  does  not 
exist  an  infinite  connected  component  a.s.  in  the  secondary  net¬ 
work  for  any  A  pt  >  0.  This  proves  part  (a). 

The  basic  idea  of  the  proof  of  part  (b)  is  to  approximate  the 
secondary  network  C?(As,  A  pt)  by  a  discrete  dependent  edge- 
percolation  model  on  the  grid.  This  discrete  dependent  edge- 
percolation  model  C  is  constructed  in  a  way  such  that  the  exis¬ 
tence  of  an  infinite  connected  component  in  C  implies  the  exis¬ 
tence  of  an  infinite  connected  component  in  G{\s,  Apt). 

Construct  the  square  lattice  C  on  U2  with  side  length  de.  Note 
that  each  site  in  C  is  virtual  and  is  not  related  to  any  node  either 
in  the  secondary  network  or  in  the  primary  network.  Next  we 
specify  the  conditions  for  an  edge  being  open  in  C,  which  is  the 
key  to  the  mapping  from  (/(As,  A  pt)  to  C. 

For  each  edge  e  in  C,  let  ( xe ,  ye)  denote  the  middle  point  of  e. 
Then  we  introduce  three  random  fields  Ae,  Be  and  Ce,  all  asso¬ 
ciated  with  the  edge  e  in  C,  where  Ce  =  AeBe  is  the  indicator 
of  the  edge  e  being  open,  Ae  represents  the  condition  (Cl)  of 
the  distance  between  two  users  for  the  existence  of  a  commu¬ 
nication  link  in  the  secondary  network  and  Be  represents  the 
condition  (C2)  of  the  spectrum  opportunity.  Specifically,  con¬ 
sider  the  Poisson  Boolean  model  B(Xs ,  rp/2,  As)  where  Xs 

9Since  a  <  rv/f  5,  every  secondary  node  in  a  square  of  is  connected  to 
those  secondary  nodes  in  the  neighboring  squares. 


E(B,e)  :=  {d(U,  B)  <  rp  —  e  for  any  infinite  connected  component  U  in  G[BC]} 


4326 


IEEE  TRANSACTIONS  ON  INFORMATION  THEORY,  VOL.  57,  NO.  7,  JULY  2011 


is  the  Poisson  point  process  generated  by  secondary  users,  then 
for  a  horizontal  edge  e,  Ae  =  1  if  the  following  two  events  (il¬ 
lustrated  in  Fig.  11)  are  true: 

(i)  there  is  an  occupied  L-R  crossing  of  the  rectangle 

\xe  -  34/4,  xe  +  34/4]  X  [ye  -  4/4,  ye  +  4/4]  in 
B(Xs,  rp/2,  As); 

(ii)  there  are  two  occupied  T-B  crossings  of  the  square  [xe  — 
3de/A ,  xe  —  4/4]  x  [ye  —  4/4,  ye  +  4/4]  and  the  square 
\xe  +  4/4,  xe  +  34/4]  x  [ye  -  4/4,  ye  +  4/4]  in 
B(Xs,  rp/2,  A s);  and  Ae  =  0  otherwise.  For  a  vertical 
edge  e,  the  definition  of  Ae  is  similar,  where  the  horizontal 
and  vertical  coordinates  are  switched  in  the  above  two 
events. 

Next,  we  define  the  random  field  Be.  For  an  edge  e  in  £, 
Be  =  1  if  Ae  =  1  and  the  following  two  events  are  true: 

(i)  there  is  no  primary  transmitter  within  distance  Rj  of  any 
secondary  node  of  the  three  crossings  in  the  definition  of 
Ae; 

(ii)  there  is  no  primary  receiver  within  distance  77  of  any  sec¬ 
ondary  node  of  the  three  crossings  in  the  definition  of  Ae; 

and  Be  =  0  otherwise.  It  follows  from  the  definition  of  com¬ 
munication  link  in  the  secondary  network  (see  Section  II-B-I) 
that  if  Ae  =  1  and  Be  —  1,  then  the  three  crossings  in 
B(Xs ,  rpj 2,  A s)  are  also  three  crossings  in  Q(X 5,  Apt). 

Let  Ce  =  AeBe ,  then  we  claim  that  the  edge  e  is  open  if 
Ce  —  1  and  e  is  closed  if  Ce  =  0.  We  observe  from  Fig.  11 
that  whether  the  edge  e  is  open  is  correlated  with  the  states  of 
the  other  edges.  This  model  C  thus  is  a  dependent  edge-perco¬ 
lation  model.  Furthermore,  as  shown  in  Fig.  12,  if  there  exists 
an  infinite  open  connected  component  in  C ,  then  those  cross¬ 
ings  associated  with  the  edges  in  the  infinite  component  in  C 
comprise  an  infinite  connected  component  in  G(Xs,  Apt).  As  a 
consequence,  by  considering  the  uniqueness  of  the  infinite  con¬ 
nected  component  in  Q (As,  Apt),  we  only  need  to  prove  the 
following  lemma  in  order  to  show  T2.1. 

Lemma  3:  Let  C(0 )  denote  the  open  connected  component 
containing  the  origin  O  in  C.  Then  given  Xs  >  Xc(rp),  3  D  > 
0,  A *PT  >  0  such  that  for  4  =  D  and  any  Apt  <  A pT,  we  have 

Pr{|(7(0)|  =  00}  >  0 

where  |(7(0)|  is  the  number  of  edges  in  C(0). 


Fig.  12.  Percolation  in  C  (thick  segments  are  open  edges  in  C  and  thin  seg¬ 
ments  are  closed  edges  in  C  and  blue  curves  are  those  crossings  associated  with 
the  open  edges). 


Proof  of  Lemma  3:  For  an  arbitrary  edge  e  in  £,  let  q  = 
Pr {Ce  =  0},  then  we  have 

q  =  Pr{(^e  =  0)  U  (Be  =  0)}  <  Pr{^4e  =  0}  +Pr{Be  =  0). 


From  the  result  on  the  crossing  probabilities  given  in  (9),  we 
know  that  when  A s  >  Xc(rp), 


Pr{A.e  =  0}  =  Pr{at  least  one  crossing  does  not  exist} 

<  I-a(Yipl)' As’  L-R 

+  [a_„((|4),As,T_B 

+  [l_„((!4),As,T_B 

— >0  as  d— >oo,Le.,  lim  Pr{Ae  =  0}  =  0. 

de  — >  oo 


Thus,  when  A s  >  A c(rp),  Ve  >  0,  3  D  >  0  such  that  Pr{Ae  = 
0}  <  |  for  4  =  D. 

Given  Ae  =  1,  let  Srt  be  the  area  of  the  region  covered  by 
the  circles  with  radii  Rj  centered  at  those  secondary  nodes  in 
the  three  crossings  and  Srj  be  the  area  of  the  region  covered  by 
the  circles  with  radii  rj  centered  at  those  secondary  nodes  in  the 
three  crossings.  Then  we  have 


Pr {Be  =  0  |  ALe  =  1} 

=  Pr{3  some  primary  transmitter  in  Srt} 

+  Pr{3  some  primary  receiver  in  Srj }. 

Since  Srt  <  +  2Rj  +  rp)  (ff  +  2Ri  +  rp)  and 

Sri  <  (^  +  2r/  +  rp)  (ff  +  2rj  +  rp),  it  follows  from 

the  basic  property  of  Poisson  point  processes  that 


Pr  {Be  =  0 
<  1  —  exp 

+  1  —  exp 


Ae  =  1} 

-Apt  +  2 Rj  + 

-Apt  (^y  +  2ri  +  (y  +  2r/  +  rJ>) 
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Obviously,  lim  Pr {Be  =  0  |  Ae  =  1}  =  0  for  fixed  de. 

A  py  — >  0 

Thus  if  we  choose  de  =  D,  then  Ve  >  0,  3  A PT  >  0  such  that 
Pr {Be  =  0  |  Ae  =  1}  <  =  for  all  A pt  <  \*PT. 

O 

It  implies  that  when  de  =  D,  for  all  Apt  <  ApT 

Pr{_Be  =  0}  =  Pr {Ae  =  0}+Pr{Be  =  0  |  Ae  =  l}Pr{Ae  =  1} 

<  =  0}  +  Pr {Be  =  0  \  Ae  =  1} 

2e 

<3- 

Thus,  for  de  =  D  and  all  Apt  <  XPT,  we  have 


A 
ft  .  : 


Primary  Transmitter 
•Secondary  User 


Fig.  13.  Counterexample  for  choosing  RPT  —  Ri.  All  the  secondary  nodes 
in  the  infinite  path  are  located  outside  those  circles  centered  at  the  primary 
transmitters  with  radii  R j ,  which  form  a  series  of  rings  surrounding  the  origin 
O  and  there  is  no  infinite  vacant  component  in  the  Poisson  Boolean  model 
B( XPTl  Rj ,  A PT)  driven  by  the  primary  transmitters. 


q  <  Pr{Ae  =  0}  +  Pr {Be  =  0}  < 


(ID 


From  Fig.  11,  we  can  see  that  if  de  >  max{4jR/  +  2rp,  4 77  + 
2 rp},  then  the  state  of  edge  e  is  only  correlated  with  its  six  ad¬ 
jacent  edges  and  it  is  independent  of  other  edges.  In  this  case, 
by  using  the  “Peierls  argument”  [25,  Chapter  1],  we  can  show 

that  if  the  probability  of  an  edge  being  closed  q  <  ^  ^  , 

then 


Fig.  14.  Infinite  path  in  the  secondary  network.  The  dashed  segments  form  an 
inner  bound  on  the  infinite  vacant  component  in  the  Poisson  Boolean  model 
driven  by  the  primary  receivers. 


Pr{|C(O)|  =  oc}>0.  (12) 

The  proof  of  the  above  statement  follows  the  one  of  [7,  Theorem 
3]  except  that  the  upper  bound  on  the  probability  of  n  edges  all 
being  closed  is  given  by  the  following  fact. 

Fact  1:  [8,  Proposition  1]:  For  any  collection  of  n 

distinct  edges  in  £,  we  have 

Pr{(Ci  =  0)  n  (C2  =  o)  n  •  •  •  n  (cn  =  o)}  <  q * 

where  C{  is  the  indicator  of  being  closed,  and  q  is  the  prob¬ 
ability  of  an  edge  being  closed. 

Thus,  by  combining  (12)  with  (11),  we  conclude  that  given 
A s  >  A c(Vp),  3  D  >  0,  A pT  >  0  such  that  for  fixed  de  = 
max{D,  4 Rj  +  2 rp,  4 rj  +  2 rp}  and  any  Apt  <  XPT 

Pr{|C(O)|  =  oo}>0. 

Notice  that  ApT  depends  on  D  which  is  chosen  according  to  the 
crossing  probability  and  is  determined  by  A5.  As  a  consequence, 
XPT  is  a  function  of  A 5,  i.e.,  XPT  =  ApT(As).  ■ 

Proof  of  T2.2:  From  the  conditions  for  the  existence  of 
a  communication  link  in  the  secondary  network  specified  in 
Section  II-B-I,  we  know  that  for  every  secondary  node  in  an  in¬ 
finite  connected  component,  there  can  exist  neither  any  primary 
transmitter  within  distance  Rj  of  it  nor  any  primary  receiver 
within  distance  77  of  it.  In  other  words,  every  secondary  node 
in  an  infinite  connected  component  must  be  located  outside 
all  the  circles  centered  at  the  primary  transmitters  and  the 
primary  receivers  with  radii  Rj  and  77,  respectively.  Thus, 
if  there  is  an  infinite  connected  component  in  the  secondary 
network,  then  an  infinite  vacant  component  must  exist  in 
the  two  Poisson  Boolean  models  B(XPt ,  Rpt ,  Apt)  and 
B(XPr,  Rpr ,  Apt)  driven  by  the  primary  transmitters  and 
the  primary  receivers,  respectively.  Here  Rpp  and  Rpr  are 


some  appropriate  radii  which  will  be  specified  later.  A  natural 
choice  for  Rpp  is  Ri ,  but  if  we  consider  the  counterexample 
given  in  Fig.  13,  then  we  can  clearly  see  that  even  if  there  is  an 
infinite  path  in  the  secondary  network,  no  infinite  vacant  com¬ 
ponent  exists  in  the  Poisson  Boolean  model  B(XpP ,  Ri ,  Apt ) 
driven  by  the  primary  transmitters.  Similarly,  counterexamples 
can  be  easily  constructed  for  choosing  Rpr  =  77. 

Suppose  there  is  an  infinite  connected  component  in  the  sec¬ 
ondary  network.  Then  we  can  find  a  sequence  of  secondary 
users  {Si,  S2,  S3, . . .}  suc^  that  they  comprise  an  infinite  path 
starting  from  Si  (see  Fig.  14). 

Assume  that  Si  and  S?;+i  (i  >  1)  are  two  adjacent  secondary 
nodes  in  the  above  infinite  path.  Notice  that  the  distance  di ?  i+ 1 
between  Si  and  S^+ 1  satisfies  di ?  i+ 1  <  rp  <  rj ,  where  the 
second  inequality  rp  <  77  follows  from  the  fact  that  the  min¬ 
imum  transmission  power  for  successful  reception  is  in  general 
higher  than  the  maximum  allowable  interference  power. 

As  we  know,  all  the  primary  receivers  must  be  outside  the 
two  circles  with  radii  rj  centered  at  S7;  and  S7;+i,  respectively, 
as  shown  in  Fig.  15.  Given  e  >  0,  consider  the  rectangle 
—  ^+1 ;  diAki  x  6]  between  Si  and  £^+1.  By  a 

simple  computation  in  geometry,  we  have  that  the  minimum 
distance  from  all  the  primary  receivers  to  the  rectangle  is 

\Aj  —  ~4  ~  e-  As  illustrated  in  Fig.  14,  it  implies  that  there 
exists  an  infinite  vacant  component  in  the  Poisson  Boolean 

model  B  ^Xpp,  \jr\  —  ^  —  6,  Apt^  driven  by  the  primary 

receivers10.  By  recalling  the  known  results  in  Section  V-A-II, 
we  thus  conclude  that  for  all  e  >  0 

Apt  <  ( 2y/ rj  -  f/4  -  e)  Ac(l). 

10This  technique  used  here  can  also  be  applied  to  the  case  when  rp  >  rIt 
where  only  the  minimum  distance  from  all  the  primary  receiver  to  the  bar  be¬ 
tween  Si  and  5i_|_i  needs  to  be  recomputed. 
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Fig.  15.  One  edge  (Si,  Si+1)  in  the  infinite  path. 


Let  6  — >  0,  then  it  yields 


XpT  - 

I  p 

The  other  term  J_  2AC(1)  in  the  upper  bound  is  obtained 
by  applying  the  same  argument  to  the  Poisson  Boolean  model 
driven  by  the  primary  transmitters. 

D.  Proof  of  Theorem  3 

Consider  the  connected  component  Ca  containing 
an  arbitrarily  chosen  secondary  user  A.  Assuming  that 
\Ca\  >  L  we  construct  a  branching  process  as  follows. 
Notice  that  if  \Ca\  >  L  then  A  must  see  an  opportunity,  i.e., 
l(A,r/,rx)  D  l(A,R/,tx)  is  true.  Call  A  the  initial  point  (or 
0-th  generation)  of  the  branching  process.  Then  the  children 
of  A  (i.e.,  the  first  generation  of  the  branching  process)  are 
secondary  users  which  satisfy  the  following  two  conditions: 

i)  it  is  within  distance  rp  of  A,  where  rp  is  the  transmission 
range  of  secondary  users; 

ii)  there  exist  neither  any  primary  receiver  within  distance  77 
of  the  secondary  user  nor  any  primary  transmitter  within 
distance  Rj  of  the  secondary  user. 

The  n-th  (n  >  2)  generation  of  the  branching  process  is  ob¬ 
tained  similarly  and  these  children  are  connected  to  their  par¬ 
ents  in  the  (n  —  l)-th  generation  of  the  branching  process  via 
bidirectional  links.  Obviously,  all  the  secondary  users  in  Ca  are 
counted  in  the  constructed  branching  process  model.  But  some 
of  them  may  probably  be  counted  more  than  once,  since  we  do 
not  exclude  the  previous  n  generations  (including  generation 
0)  when  we  consider  the  n-th  generation.  Thus,  this  branching 
process  gives  us  an  upper  bound  on  the  number  of  secondary 
users  in  Ca-  It  follows  that  if  the  branching  process  does  not 
grow  to  infinity  wpp.,  then  there  does  not  exist  an  infinite  con¬ 
nected  component  a.s.  in  G(\s,  ^pt)>  due  to  the  stationarity 
of  the  heterogeneous  network  model.  Since  the  conditional  av¬ 
erage  degree  p  is  the  average  number  of  offspring  for  every  gen¬ 
eration,  the  necessary  condition  follows  immediately  from  the 
classic  theorem  for  branching  processes11  [1,  Theorem  2.1.1]. 

nThe  lower  bound  on  the  conditional  average  degree  can  be  improved  by 
constructing  a  multi-type  branching  process.  Specifically,  when  we  consider  the 
children  of  a  secondary  user,  we  categorize  them  into  different  types  based  on 
their  distances  from  the  secondary  user.  A  similar  construction  has  been  used 
to  obtain  the  lower  bound  on  the  critical  density  of  Poisson  Boolean  models  [1, 
Sec.  3.9].  This  multi-type  branching  process  is,  however,  difficult  to  analyze  and 
requires  substantial  further  investigation. 


E.  Proof  of  Theorem  4 

From  the  construction  of  the  dependent  site-percolation 
model  C ,  we  know  that  the  existence  of  an  infinite  occupied 
component  in  C  implies  the  existence  of  an  infinite  connected 
component  in  G(\s,  Apt)-  Then  in  order  to  obtain  a  sufficient 
condition  for  the  connectivity  of  the  secondary  network,  it 
suffices  to  find  a  sufficient  condition  for  the  existence  of  an 
infinite  occupied  component  in  C. 

Let  p  be  the  probability  that  one  site  is  occupied.  Then  based 
on  the  definition  of  the  upper  critical  probability  pc  of  C,  we 
have  that  if  p  >  pc,  an  infinite  occupied  component  containing 
the  origin  exists  in  C  wpp.  It  implies  that  if  p  >  pc,  there  exists 
an  infinite  connected  component  in  the  secondary  network  wpp. 
Since  the  event  that  there  exists  an  infinite  connected  component 
in  the  secondary  network  is  invariant  under  the  group  of  shift 
transformations,  it  follows  from  the  ergodicity  of  the  heteroge¬ 
neous  network  model  (see  Lemma  1)  that  if  p  >  pc,  there  exists 
an  infinite  connected  component  in  the  secondary  network  a.s. 

Let  Fk(Bo)  denote  the  event  that  there  exist  exactly  k  sec¬ 
ondary  users  in  the  associated  dashed  square  Bo  of  a  site  O  and 
let  Gi(Bo)  (1  <  i  <  k)  denote  the  event  that  the  secondary 
user  Ai  sees  an  opportunity.  Then  based  on  the  definition  of  oc¬ 
cupied  site  in  C,  we  obtain  a  lower  bound  on  p  as  follows: 


p=  ^Pr{Ffe(B0)}Pr{ui1Gi(i?o)} 

k=  1 

OO 

>  ^Pr{Ffe(Bo)}Pr{Gi(So)} 

k= 1 


=  Pr{G1(Bo)}^Pr{Ffc(Bo)} 

k= 1 


=  [1  -  exp(— Asc^)]Pr  jl  (  4i .  /  /.  rx )  n  1 


1  — exp 


exp{-ApT7T  [. Rj+rj-I(Ri ,  Rp,  r/J}, 


In  the  last  step,  Pr  j  I  (A\ ,  77,  rx)  fl  0  (A\ ,  i?/,  tx)  j  has  been 
obtained  by  setting  the  distance  d  =  0  in  the  expression  for  the 
probability  of  an  unidirectional  opportunity  between  two  sec¬ 
ondary  users  with  distance  d  apart  given  in  [16,  Proposition  1]. 

VI.  Conclusion  and  Future  Directions 

We  have  studied  the  connectivity  of  a  large-scale  ad  hoc  het¬ 
erogeneous  wireless  network  in  terms  of  the  occurrence  of  the 
percolation  phenomenon.  We  have  introduced  the  concept  of 
connectivity  region  to  specify  the  dependency  of  connectivity 
on  the  density  of  the  secondary  users  and  the  traffic  load  of 
the  primary  users.  We  have  shown  several  basic  properties  of 
the  connectivity  region:  the  contiguity  of  the  region,  the  mono¬ 
tonicity  of  the  boundary  and  the  uniqueness  of  the  infinite  con¬ 
nected  component.  We  have  analytically  characterized  the  crit¬ 
ical  density  of  the  secondary  users  and  the  critical  density  of 
the  primary  transmitters;  they  jointly  specify  the  profile  of  the 
connectivity  region.  We  have  also  established  a  necessary  and  a 
sufficient  condition  for  connectivity,  which  give  an  outer  and  an 
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(C) 


Fig.  16.  Illustration  of  SI(tJrllr2)  (the  common  area  of  two  circles  with  radii 
r  i  and  r2  and  centers  t  apart),  SU2(t,r1,r2)  (the  union  area  of  two  circles  with 
radii  r1  and  r2  and  cents  t  apart)  and  SI2(r707  Rp,t,rj)  (the  intersection  area 
between  one  circle  with  radius  Rp  and  the  union  of  the  two  identical  circles 
with  radii  rT). 


Proposition  1:  Let  A s  and  A pt  be  the  density  of  secondary 
users  and  primary  transmitters.  Let  77  and  Rj  be  the  interfer¬ 
ence  range  of  the  secondary  and  primary  users  and  rv  and  Rp 
the  transmission  range  of  the  secondary  and  primary  users.  Let 
rp,  77,  Rv.  Rj)  denote  the  probability  of  a  bidirectional 
opportunity  between  A  and  B  conditioned  on  the  event  that  A 
sees  an  opportunity.  Then  the  conditional  average  degree  g  of  a 
secondary  user  is  given  by 


(Astt r-p)  •  5(ApT,rp,rj,i?p,i?/) 
rrp  2 1  \  r 

=  AsTiTp  /  ^exp  <  -  Apt  7r(rf  +  Rj  +  I(Ri,  Rp,  rj)) 
Jo  rp  [  1 

-  Sj(£,  77,  77)  -  £/(£,  jRj,  i?/) 


inner  bound,  respectively,  on  the  connectivity  region.  Further¬ 
more,  by  examining  the  impacts  of  the  secondary  users’  trans¬ 
mission  power  on  the  connectivity  region  and  on  the  conditional 
average  degree  of  a  secondary  user,  we  have  demonstrated  the 
tradeoff  between  proximity  and  spectrum  opportunity.  In  es¬ 
tablishing  these  results,  we  have  used  techniques  and  theories 
in  continuum  percolation,  including  the  coupling  argument,  er- 
godic  theory,  the  discretization  technique,  and  an  approximation 
using  a  branching  process. 

To  highlight  the  interaction  between  the  primary  users  and  the 
secondary  users  in  a  heterogeneous  network,  we  have  ignored 
the  fading  effect  and  the  mutual  interference  between  secondary 
users.  If  we  take  into  account  these  factors,  then  the  received 
signal  to  interference-plus-noise  ratios  at  two  secondary  users 
will  replace  the  distance  between  them  in  condition  (Cl)  for  the 
existence  of  a  communication  link  between  them.  This  will  re¬ 
sult  in  a  random  connection  model  (RCM)  with  correlated  links, 
where  the  correlation  between  links  is  due  to  the  mutual  inter¬ 
ference  and  condition  (C2)  on  the  presence  of  the  bidirectional 
opportunity.  There  will  be  a  tradeoff  between  proximity  and  mu¬ 
tual  interference  in  addition  to  the  tradeoff  between  proximity 
and  opportunity  and  the  combination  of  these  two  tradeoffs  will 
complicate  the  characterization  of  the  connectivity  of  the  sec¬ 
ondary  network.  But  since  the  RCM  shares  several  basic  prop¬ 
erties  (e.g.,  ergodicity  and  existence  of  the  critical  density)  with 
the  Boolean  model  used  in  this  paper,  we  expect  that  most  of  the 
results  established  here  (e.g.,  connectivity  region,  critical  den¬ 
sities)  can  be  extended  to  the  RCM. 

Appendix  A 

Expression  for  Conditional  Average  Degree 


Sj(t,  77,  T2)  the  common  area  of  two  circles  with  radii  77  and 
77  and  centers  t  apart  [see  Fig.  16(a)]  and  Sjj2(t ,  n,  77)  is  the 
union  of  two  circles  with  radii  77  and  77  and  centers  t  apart  [see 
Fig.  16(b)].  Si2(r ,  0,  Rp,  t,  77)  is  the  intersection  area  between 
one  circle  with  radius  Rp  and  the  union  of  the  two  circles  with 
both  radii  77  [see  Fig.  16(c)].  For  Sj2(r,8,Rpti9ri),  the  two 
identical  circles  are  centered  t  apart  and  the  other  circle  is  cen¬ 
tered  at  (r,  6),  where  the  middle  point  of  the  centers  of  the  two 
identical  circles  is  chosen  to  be  the  origin  O. 

Expressions  for  I(Rj,  Rp,rj)  and  Sj(t,  77, 77)  can  be  ob¬ 
tained  in  explicit  form,  which  can  be  found  in  [1,  Appendix  A]. 
The  expression  for  Sj2 (r,  <9,  Rp ,  t,  77)  depends  on  the  expres¬ 
sion  for  the  common  area  of  three  circles  which  is  tedious  and 
is  given  in  [26].  By  applying  the  basic  property  of  the  exponen¬ 
tial  function  to  (Al),  we  can  easily  show  that  g(-)  is  a  strictly 
decreasing  function  of  A  pt- 

Proof:  Let  Fk  (A)  denote  the  event  that  there  exist  exactly  k 
neighbors  of  a  secondary  user  A.  We  thus  have  the  first  equation 
shown  at  the  bottom  of  the  next  page. 

When  k  —  0,  it  is  obvious  that  deg(A )  =  0.  When  k  >  0,  let 
Bi  be  a  neighbor  of  A  and  1  pi  an  indicator  function  for  Bi  such 
that  1  Bi  =  1  if  D (Bi,  rj ,  rx)  fl  I ( B, ,  Rj ,  tx)  occurs  and  1  Bi  — 
0  otherwise.  Then  by  considering  the  statistical  independence 
and  equivalence  of  the  k  secondary  users,  we  have  the  second 
equation  shown  at  the  bottom  of  the  next  page.  It  follows  that 
we  get  (A2),  as  shown  at  the  bottom  of  the  next  page. 

According  to  the  definition  of  spectrum  opportunity, 
Pr{D(A,  77,  rx)  fl  l(A,  jRj, tx)}  can  be  obtained  by  setting 
the  distance  d  —  0  in  the  expression  for  the  probability  of  a 
unidirectional  opportunity  between  two  secondary  users  with 
distance  d  apart  given  in  Proposition  1  in  [16]. 


An  expression  for  the  conditional  average  degree  g  of  a  sec¬ 
ondary  user  is  presented  in  the  following  proposition. 


Pr{I(A,  77,  rx)  fl  I  (A,  Rj ,  tx)} 

=  exp[-XPTTT(rj  +  Rj-  I{RI,Rp,n ))]  (A3) 
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Next,  we  derive  the  expression  for  the  probability 
of  a  bidirectional  opportunity,  i.e.,  Pr{l(Si,  77,  rx)  D 
I ( /?  1 ,  ,  tx) nI(A,r/,  rx)  fl  I  (A,  Sj,  tx)},  which  depends 
on  the  location  of  B\  only  through  its  distance  to  A.  Since  Si 
is  uniformly  distributed  within  distance  rp  of  A,  the  density 
function  of  the  distance  t  between  B\  and  A  is  given  by  ^ 

rp 

for  0  <  t  <  rp.  In  this  case,  the  probability  of  a  bidirectional 
opportunity  can  be  written  as  shown  in  (A4),  also  shown  at 
the  bottom  of  the  page,  where  the  integrand  can  be  written  as 
shown  in  (A5),  at  the  bottom  of  the  page. 


Next,  we  compute  the  two  probabilities  in  (A5)  one  by  one. 
Since  the  primary  receivers  admit  a  Poisson  point  process  with 
density  \pp,  we  have 


Pr{D(Si,r/,rx)  n  I(A,r/,rx)|  d(Bi,A)  =  t} 

=  exp[-ApT(27rrf  -  Sj(t,  rj,  r/))]  (A6) 

where  Sj(t ,  77, 77)  is  the  common  area  of  two  circles  with  both 
radii  77  and  centers  t  apart  [see  Fig.  16(a)]. 


/x  =  E  [deg(A)  |  I  (A,  77,  rx)  fl  D  (A,  Rp  tx)] 

=  Ek[E[deg(A )  |  l(JL,.Tj,rx)  n  I(A,i?/,tx)  n  Fk(A)]] 

00  ( \  2\ ^ 

=  y^e_As7rrp — — E[deg(A)  |  l(A,rj,rx)  n  D(A,i?/,tx)  n  Fk(A)\ 

k= 0 


E  [deg(A)  |  I  (A,  77,  rx)  fl  I  (A,  Sj,  tx)  fl  Fk(A)\ 

k 

=  ^  E[1  Bi  |  KA rj, rx)  n  1(1, i?/,  tx)] 

i=l 

=  kE[lBi  |  B(A,r/,rx)  fl  D(A,  S/,tx)] 

=  &Pr{l(Si,  rp  rx)  fl  l(Si,S/,tx)  |  l(A,r/,rx)  fl  D (A,  S/,tx)} 

^  Pr{l(Si,  77,  rx)  fl  I  (Si,  Sj,  tx)  D  l(A,  77,  rx)  fl  D(A,  Sj,  tx)} 
Pr{l(A,  77,  rx)  fl  l(A,  Sj,  tx)} 


g  =  XsTTTp 


Pr {I (Si, rj, rx)  fl  l(Si,S/,tx)  fl  l(A,r/,rx)  fl  l(A, S/,tx)} 
Pr{l(A,  77,  rx)  fl  I  (A,  Sj,  tx)} 


(A2) 


Pr{l(Si, 77, rx)  fl  I(Si,S/,tx)  fl  l(A,r/,rx)  fl  l(A,  S/,tx)} 


Jo 


p  2 1 


-^-Pr{D(Si,  77,  rx)  fl  I  (Si,  Sj,  tx)  fl  I  (A,  77,  rx)  fl  l(A,  Sj,  tx)|  d(Si,  A)  =  t}dt 


(A4) 


Pr{l(Si,  77,  rx)  fl  I  (Si,  Sj,  tx)  fl  D  (A,  77,  rx)  fl  l(A,  Sj,  tx)  |  d(Si,  A)  =  t} 

=  Pr{I(Si,  Sj,  tx)  fl  I  (A,  S/,  tx)  |  I  (Si,  77,  rx)  fl  I  (A,  77,  rx)  fl  d(Si,  A)  =  £} 
•  Pr{l(Si, 77,  rx)  fl  D(A,  77,  rx)  |  d(BpA)  —  t} 


(A5) 
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Pr{I(i?i,  Ri,  tx)  fl  I  (A}  Rp  tx)  |  D  ( /?  i ,  rj.  rx)  fl  I  (A,  77,  rx)  D  d(BpA)  =  t} 

Si2(r,6,Rp,t,ri) 


=  exp  ^  -Apt  J  J 

Su2(t,Ri  ,Ri) 


1  - 


irRj 


-rdrdO 


=  exp  <  -XPT 


2pR2j  -  Si(t,  Ri ,  iZ/)  -  J  J 


Si2(r,0,Rp,t,ri) 


Su2(t,Ri,Ri) 


TtRp 


rdrd# 


(A7) 


/i  =  A^TTr^ 


2b 

^exp 

rP 


-Xpt  [  vrr?  —  2 r?  arccos  f  — )  +  t 
\  V>i/ 


d  t. 


rrv  2 1  ( 

M  3  exp  ^  -APT 

Jo 


■xrZi-2rj[--  —  )+t\lrj-f- 


frp  2t 

<  XS7rr2  /  exp(-APT£rj)d£ 

Jo  rp 

(  2  2  2/3 

=  Astt  P  2  -  -rg — 2  exp(-APT/3r|)  -  —t—  exp(-APT/?rf ) 

\^PTrI  ^PTr  I  App 


2\stt  ,  .  2 

< 

App 


d  t 


Pr{l(i3i,  77,  rx)  D  iZj,  tx)  fl  I(A,  77,rx)  fl  I  (A, i2/,tx)  |  d(i?i,  A)  =  £} 
=  exp  |  -  Apt  [2x(t?  +  -  Sj(t, rr,rr)  -  Si(t, Rr,  Rj ) 


~  If 

Su2(t,Ri  ,Ri) 


Sj^pOPppri) 

7 rR2p 


rdrdO 


(A8) 


Let  Xpp  denote  the  Poisson  point  process  formed  by 
primary  transmitters.  If  we  remove  from  Xpp  primary  trans¬ 
mitters  whose  receivers  are  within  distance  77  of  Bi  or  A, 
then  it  follows  from  the  Coloring  Theorem  [17,  Chapter  5]  that 
all  the  remaining  primary  transmitters  form  another  Poisson 


Si2('r,6,Rp,t,ri) 

7 xRl 


where 


point  process  with  density  Xpp 
S12 (r,  0 ,  Rp,  f,  rj)  is  the  area  of  the  circle  with  radius  Rp  and 
centered  at  (r,  0)  intersecting  the  two  circles  both  with  radii 
77  and  centers^  apart  [see  Fig.  16(c)].  We  thus  have  (A7),  at 
the  top  of  the  page,  where  Rp  Ri)  is  the  union  of  two 

circles  with  both  radii  Rj  and  centers^  apart  [see  Fig.  16(b)]. 

Substitute  (A6),  (A7)  into  (A5),  we  have  (A8),  as  shown  at 
the  top  of  the  page.  The  expression  for  the  conditional  average 
degree  /i  thus  follows  by  plugging  (A8)  into  (A4)  and  then  (A3), 
(A4)  into  (A2). 


Appendix  B 
Proof  of  Corollary  3 

From  [1,  Appendix  A]  and  Figs.  16(b)  and  (c),  we  know  that 
when  n  >  Rp  +  Ri 


I(RI,Rp,n)=Rj  (Bl) 

//  e=SB2(t,R„R,) 

Su2(t,Ri  ,Ri) 

=  2pR2j  —  Sj(t ,  iZj,  Rj ). 

(B2) 


Substitute  (Bl)  and  (B2)  into  (Al),  we  have 
irp  2 1 

H  =  \sir r\  /  -j  exp[— ApT(7rrf  —  S/(L  rj,  r/))]d£.  (B3) 
3o  rp 
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Plugging  the  expression  for  Sj(t,  77.  77)  [1,  Appendix  A]  into 
(B3)  yields  the  equation  shown  at  the  top  of  the  previous  page. 

By  applying  the  inequality  arccos(rr)  <  xforO  <  x  <  1, 
we  have  the  equation  at  the  top  of  the  previous  page,  where  we 
have  assumed  that  rp  =  [3rj  (0  <  (3  <  1)  under  the  disk  signal 
propagation  and  interference  model.  Since  77  oc  (ptx)1^,  we 
arrive  at  Corollary  3. 
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